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COMPUTER PROGRAM FOR CALCULATING FLOW DISTRIBUTION 
IN A RADIAL-INFLOW TURBINE 

by Theodore Katsanis 


SUMMARY 

A FORTRAN computer program for flow analysis of a radial inflow gas 
turbine is given. The program obtains a meridional solution on the mean 
surface between the blades, followed by solutions on hub, mean, and shroud 
blade-to-blade surfaces, in- a single computer run. Suggestions for modify- 
ing the program for use with other types of turbomachines are given. Tech- 
niques for overcoming convergence problems are discussed. 


INTRODUCTION 

A method of flew analysis for any turbomachine is summarized in refer- 
ence 1, with the details being presented in references 2 and 3. Reference 2 
gives the method applied to the meridional plane. Reference 3 gives the 
method for obtaining a blade-to-blade solution, using information obtained 
from the meridional plane solution. The FORTRAN program for each solution 
as applied to a radial inflow gas turbine is given in the appropriate ref- 
erence. 

As stated in reference 3, a FORTRAN program has been written which com- 
bines both programs. The combined program obtains first the meridional 
solution, followed by three blade-to-blade solutions (hub, mean, and shroud), 
in a single computer run. Either the meridional solution or blade-to-blade 
solution may be obtained separately. 

This report presents this FORTRAN program, with sample output for the 
example rotor of reference 1. Also suggestions are given for modifying the 
program' for other types of turbomachines, and for overcoming convergence 
problems which may be encountered. 


PROGRAM VARIABLES 

The variable names used in the combined program are the same as listed 
in reference 2 and 3 with the addition of the following variables. 

THH(I) 6 coordinate of mean blade surface at quasi -orthogonal 

along hub 

THM(I) 0 coordinate of mean blade surface at I quasi -orthogonal 

along mean streamline 
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THS(-I) 

9 

coordinate of mean blade surface at I^ quits i -orthogonal 
along shroud 

PION 

tf/N 

TPION 

2*/N 

THHI(l) 

9 

tb 

coordinate of suction surface of complete blade at I 
quasi -orthogonal along hub 

THMT(I) 

0 

coordinate of srction surface of complete blade at I** 
quasi -orthongonal along mean blade -to-b lade stream surface 

THSI(I) 

9 

coordinate of suction surface, of complete blade at I^* 1 
quasi -orthogonal along shroud 

THHKH(I) 

9 

*f*Vi 

coordinate of pressure surface of splitter blade at I 
quasi -orthogonal along hub 

thmech(i) 

0 

til 

coordinate of pressure surface of splitter blade at i 
quasi -orthogonal along mean blade -to-b lade stream surface 

THSKH'I) 

0 

coordinate of pressure surface of splitter blade at I th 
quasi -orthogonal along shroud 

THHKP(l) 

0 

coordinate of suction surface of splitter blade at Tth 
quasi -orthogonal along hub 

THMKP(l) 

9 

coordinate of suction surface of -splitter blade at I** 1 
quasi -orthogonal along mean blade-to-blade stream surface 

THSKP(.I) 

0 

coordinate of .suction surface of splitter blade at I** 
quasi -orthogonal along shroud 

THHKMX(I) 

9 

coordinate of pressure surface of complete blade at I th 
quasi -orthogonal along hub 

THMKMX(I) 

0 

coordinate of pressure surface of complete blade at I th 
quasi -orthogonal along mean blade -to -blade stream surface 

THSKMX(l) 

0 

coordinate of pressure surfsce of complete blade at I^* 1 
quasi -orthogonal along shroud 


MODIFICATIONS REQUIRED FOR OTHER APPLICATIONS 

The program as written can "be used for a radial inflow gas turbine. 
It should be emphasized that the listing of the combined program is being 
published as a guide to the technique used in programming the quasi- 
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orthogonal method of f ? analysis. If it is desired to obtain a program 
for analysing any type of turbomachine other than a radial inflow gas tur- 
bine, several modifications as indicated below will have to be made to the 
program. This will require a reasonable understanding of the program and 
the equations involved. With this understanding, the required modifications 
can be easily made. Modifications which are necessary for some particular 
cases are listed below: 

(1) Changes required for a compressor are: 

(a) dotation of coordinates for use in calculating spline fit 
curve for the meridional streamlines. It is assumed that the inlet 
is to the left with the positive direction to the right, and that 
the quasi -orthogonals are numbered starting with 1 at the inlet. 

The coordinates for the streamlines must be rotated 45° in the di- 
rection opposite to that in the present program. The modifications 
needed sire to statements for calculating AB(l' and AC(l) near state- 
ments 160 and 390, and for calculating AL(l,iC) four statements after 
statement 160. - 

(b) The sign of CEP as calculated just before statement 150 must 
be reversed. 

(2) Changes required for a pump or liquid turbine. 

Statements involving temperature or density must be modified 
to eliminate temperature, and to allow for constant density. 

(3) Changes required for axial flow machines. 

Statements for rotating streamline coordinates for use in calculat- 
ing the spline fit curve for the meridional streamlines should be 
modified to eliminate this rotation. The statements involved are 
the same as those mentioned under la above. 


CONVERGENCE 

It may be found that with some geometries there is a convergence prob- 
lem. This is helped by using as few and as evenly spaced quasi -orthogonals 
as possible, consistent with a reasonable accuracy. The number of stream- 
lines used does not appreciably affect the. number of iterations. Of course, 
it does affect the computer time required per iteration. Another factor 
which helps convergence is to reduce the streamline correction factor, CORFAC, 
(see ref. (2), p. 23 and fig. 5). Of course, this reduces the rate of con- 
vergence, but may result in convergence when a solution could not be obtained 
otherwise. 

However, in some cases when the initial streamline estimate is poor it 
may be necessary to make the streamline correction factor, CORFAC, so small 
that the convergence rate is too slow. One way of improving the initial 



streamline estimate is to use a large streamline correction for the first 
iteration, say_ CORFAC = 1* This will bring the streamlines fairly close to 
the correct position, but will result in rather uneven streamlines. If the 
streamlines are then smoothed sufficiently, the convergence will be generally 
much better than when working from the original position. It has not been 
found helpful to make this large streamline correction more than once. There 
are several mathematical techniques for smoothing (ref. 4). 

Another problem that may be encountered with compressible fluids is that 
calculations may indicate choking weight flow less than the desired weight 
flow in the early interations. This usually leads to a problem in conver- 
gence. However, this may be overcome by obtaining a solution based on a 
slightly lower weight flow, followed by a solution based on the desired weight 
flow. Possibly more than one increment of weight flow increase will be re- 
quired to increase the flow up to the desired weight, or perhaps the desired 
weight flow will prove to be actually in excess of the choking weight flow. 
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PROGRAM LISTING 


$I8FTC FIXED OECK 

SUBROUTINE F I.NEOI HUB » MEAN , STROUD) 

CALCULATION 3F VELOCITY AMD PRESSURE DISTRIBUTION IN A RADIAL FLO* 
TUR83 MACHINE 


COMMON S8W,KMX ,MX,WT,XN ,3AM, AR, TYPE, BC3P, TEMP,ALM,RHO,TOL£R, 
to 1 PLOSS, NPRT, I TER1 »3ETIN* WTOLER, THHl ,THHKH,THHKP,THHKMX,THM1, 

to 2 THMKH,'"HMKP»TH«KMX,THS1 ,THSKH t 'TSKP,TTSKNX,Z ,B ,AB,AC,AO, 

J, 3 RUNO, M)'BL »P IQN »W,8ETAI »WTR,CURV 

DIMENSION AL(21,21),9ETA(21,21) ,CAL(2I,2l» ,CBF TAI21 ,21 » , 

1CURVI 21, 22), DM (21,211, PR S(2l ,21) ,R(2) ,215 , Z (2 1 ,21 ) ,SM(21,21 ) , 

2S A( 21, 21), SB (21, 21), SC (21, 21?, SD (21 ,21) ,SAL(2i ,21) ,SBETAI21,21 f, 
3TN(21,21),TT(21,21),WAI21,21) ,WTR(21 ,21) 

OIMEISION AB( 22),AC(22),AO(22),BA(21) ,DELBTA(21) ,ORDH(2l>» 

IDTDRl 21) sOTDZ(2l) ,3WMDM(21) ,DWTOM(2H ,RH (2 1 ) ,RS ( 211 ,7H(21 ),ZS( 21 ) , 
2THTA(2i),HTr L (21 ) ,XR C21 ) , XT« 21 ) ,XZ (21 ) ,BETAI (3) ,AA(3) 

DIMENSION THH( 21) ,THM(21 ) , THS(21 ) ,THH1 (21) ,THHKH(21 ) ,THHKP( 21) , 

1 THHXNX( 2l),THMl(?i) »THMKH 5 21) ,THMKP(21) ,TTMKMX(21 K, THS1 ( 2 1 ) . 

2 THS< H( 21 ) , THSKP ( 21 ) , THSKMX (21 ) 

INTEGER RUNO, TYPE ,BCDP, SRW,HUB , SHROUD 
RUNOs-O 

10 READ ( 5, 1010)MX,KMX,MR,MZ,W,HT,XN,GAM,AR 

ITNO a 1 

RUN0=RUN3+1 

WRITE (6, 1020) RUNO 

WRITE (6,1010 )MX»KM: ,MR,MZ,W,WT,XN,GAM,AR 

READ (5, 1010) TYPE, 3C0P,SRW,MXBL, TEMP, ALM,RHO, TOLER, PLOSS,WTOLER 

WRITE (6, 1 010 ) TYPF.3CDP, SRW,MXBl, TEMP, ALM,RHO, TOLER, PLOSS,WTOLER 
READ (<S, 1010)MTHTA ,NPRT,I TER ,NULL ,SFACT,ZSPLI T., BET IN, CORF AC 

WRITE! 6, 1010)MTHTA,NPRT,ITER,NULL,SFACT,ZSPLIT,BETIN,CORFAC 
ITER1 = ITE1 

READ( 5, 1030) (ZSJ I), 1*1, MX) 

WRITE(6, 1030)(ZS( t ),I*1,MX» 

READ( S, 1030)(ZH( I ), 1*1, MX) 

WR I TFi 6, 1030 )(ZH(T ),I*1,MX) 

R EAD( 5, 1030) ( RS( I ), 1*1, MX) 

WRITF(6, 1030)(RS(I ),I=l,MX) 

READ! 5, 10301 ( RH( I ) , 1*1 ,MX) 

WRITE(6, 1030)(RH( I ),T=l,MX) 

DO 20 1*1, MX 

, zs( n*zs( n/ 12 . 

ZH( I ) *7H( I ) / 12. 

RSI I)*RSl 1 5/12. 

20 RH( I)*RH( I)/12. 

• IF( TYPE.NE *0 ) GO TO AO 

WA( 1,1) » WT/RHO/(ZSIl)-ZH(l) ) /3. 1 A / (RSI 1 ) +RH( l ) ) 

DO 30 1*1, MX 

ON ( I, KMX ) *SORT( ( ZS( I)-ZH(I) )**2* (RS ( 1 ) -RH( I ) ) **2) 

DO 30 K*1»KMX 

ONI I,K)«FL0AT(K-1)/FL0AT(KMX-1)*DN(I ,KMX) 

WA(I,K)*WA(1,1) 

Z( l,K)*DN( I,K)/DN(I,KMX)*(ZSm-ZH(I)M>ZHII) 
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30 R( It < )=DN( It KI/ 3 N( I,KMX)*(RS(I ) — R-l ( T > ) MH( II 
GO TO 50 

. AO IF( TYPE .ME . 1 ) GO TO MS 

CALL BCR EA[. ( ON ( 1 , l > , DN ( 21 , 2 1 >> 

CALL 6CREA0 I WA( 1, 1 I , WA ( 21 *211 1 
CALL BCRPAO ( Z ( l , 1) , Z 1 21 , 211 1 
CALL BCREAO (Rf 1, M ,R1 21 ,21)} 

WRITF (6,1040) - 

50 READ ( 5, 1030H THTA (I ) , I *1 ,MTHTA ) 

•(RITE ( 6, 1030) ( THTA ( I ) ,1 *1 ,MT-ITA> 

READ (5, 1030MXT(IM-1 ; NT4TAI '' 

WRITE! 6, 1030 MXT( I ),l=l,MT-(TA) 

00 60 K*1,MR 

READ (5,1030J(TN(I,K>,I*1»MZ> ' ~- 

60 WRITE (6,1030HTN(I,K'/,I»1,MZ) - ■- 

READ (5,1030HXZ(I),1»1,MZ) 

WRITE (6, 1 030) ( XZ( I ),I*1,MZ) 

READ I5,1030HXR(n,f*.l,MRl 
WRITE (6, 1030KXRU 1,1 = 1, HR) 

END OF IN n UT STATEMENTS 

Sr ' AL rr!M«T^Te E 1NCHES Til FfcCT ANr ' PSr T0 IH/SQ ET , INITIALIZE, CALCJLAT 


70 00 90 <*I,MR 
DO 30 1*1, MZ 

80 TN{ I,K I « TN(I,K>/1’« 
90 XR( K 1 * XR(<)/12. 

DO >00 I "1 ,MZ 

ioo xzm * xz( i i/i2. 

00 110 K»1,KMX 

no SNU»KI*0. 

BA(l)»0. 


DO 120 K«2,KMX 

120 BAlK) « FL0AT(R-1>*WT/FL0 ATvX« -1 > 
00 130 I»1,MX 
130 DNl-I, U»0, 

00 140 I«1,MTHTA 

140 xrm-xTd)/i2. 

ROOT - SORT! ’.O) 

145 CONTINUE 

TOLER •TOLER/ 12. 

Z SPLIT • ZSPLIT/12. 

PL OSS -PL OS $*144. 

Cl ■ SQRT( GAM* AR* TEMP > 


I Tt ! 6,1050 1 Cl 
KWXM1 ■ KMX- l 
CP»AR*GAM/(GaM-U) 

EXPOV « l./IGAM-l.) 

BETIN ■ BET IN/57, 2,9577 
ZINLFT • (ZS(l)*ZH(l))/2. 

RINLET » (RS(lHRH(in/ 2 . 

CALL LIMINT (ZIN16T,RINLET,X2,XR,TN,21 ,21 ,T) 

RB • RINLfT*EXP(-. 71*(2. *3,141 59/IXN*SFACT»-T/RINLE1 
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WRITE (6,1030) R9 

CEF = -SIN(BE»"IN)/CO$(BETIN)/RINLET/(RINLET-RB) **2 
ERR0R=100000* 

BEGINNING OF LOOP FOR ITERATIONS 


150 IF< ITEft.E-%0) WRITE (6,1060) ITNO 
IFCITEx.EQ.O) WRITE (6,1070) 

ERROR 1=ERR0R 
ERR0R*0. 

START CALCULATION OF PARAMETERS 


to 

I 

W 


00 230 K = 1 , KM X 
DO 160 1=1, MX 

ABC 1 1 * (Z(I,K)-R( I ,K ) ) /ROOT 

160 ACm«(4U».)+R<I(Kn/R00T " * . 

CALL SPLINE (A8,AC,MX,AL(1,K) ,CURV ( 1 , K) ) 

00 170 1= l, MX 

CURVt I,K)=CURV( ? tK)/( l.+AL(I,K) **2) **1 « 5 
AL( I, K ) « ATAN( AL( ItK) )-* 78M398 
CAL ( If K 1 = C0$( AL ( !«K) ) 

170 SAL ( I ,K ) = S IN ( AL ( I, K ) ) 

00 180 1=2, MX i 

180 SMC I, K ) = $M( I- 1,K ) + SQRT( (Z(I ,10-2(1-1 ,K) R(IfK)-RII-ifK)#** 

1 2) 

190 CALL SPLDERC XT( 1 ) , THTA ( 1 1 , MTHTA ,Z (1 ,K) , MX, OTOZ( 1 ) ) 

DO 220 1=1, MX 
T = 0. 

IF( I.LEfMXBL) CALL L ININTC Z( I ,K) ,RU ,K) ,XZ ,XR,TN, 21 , 21 ,T ) 

IF(R( I,K) .lF.RB)G3 TO 200 
OTDR< I )»CFF*( R ( I,K)~R8)**2 
GO TO 210 
200 OTDRUMO. 

210 TO=R( I f O*DTDR<I) 

TP = Rl !fKI*DTDZ(II 
TT( I,K)=T*SQRT(1.+T2*TP) 

BET At 1,0 =ATAN< TP*C4L ( I * K) +TQ*SAL( I »K) ) 

S8ETA( T , K ) - S INC BETA ( I ,K M 
CBFTAdfO C3$( BETA ( I ,K) ) 

SA( I, K)=C8ETA( I ,K V**2*CAL 1 1 f K)*CURV (I f K) -SBb T A( I , K) W/RII t KM- 
1SAK If<)*CBcTA(IfK1*S3FTA(I f K)*OTTR(I) 

SCC l f K)=-SAL( If<H*C8FTA (I , K ) **2*CURV (I tKMSAl ( T # K) *CRET ACT »K> 
i*SRET4( I,<)*DTDZ(I) 

AB( n=WA( IfK )*CBETAU,K) 

220 ACUMWAUfKMSBETAUtM 

CALL Sr>l!NFC SM ( 1 ,K ) f AB , MX, OWMOIM ,A 0) 

CALL SPL IN 6 C SM(t ,K 5 ,AC ,MX , OWTOM.AO! 

IFC ( IT£R.LF*0).ANP*(M0D(K“1 , NPRT)«EQ*0) ) WRITE (6,1080) K 
DO 230 1=1, MX 

SB t 1 , K )=$AL ( I»K)*CBETA(I ,K) *f)WMOM (I)-2#^W*SBETA(I f K MOTDRf I )* 

1R( I,< MC3ETAC t»K)*(DWTOMU M2.*W*SAL(I , K) ) 

S0( I, KMCALt I,K)*r.BFTA(tfK)*3WMDMU MOTOZUI* 

IRf I,K )*CBFTA( I ,K )♦( OWTOMII M?.*W*SAL (I , K) ) 

I F( ( I TER *GT #0 ) «TR * C MOD( K- 1 • NPRT) *N£ # 0) ) 30 TO 230 
A* AL ( ItK )* *>7*29577 
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B= SWII,KI*12. 

E= TTII,K>*12. 

G=6ETAI I,X 9*57.29577 

WRITF (6,10<»n| A,CU?VU ,K) ,B,G,F, S Af I } K‘. , SB ( I , K» , SC ( I ,K 9 , SOI I, K ) 
230 CONTINUE 
C 

C END OF LOOP - ®ARANETER CALCUtATION 

C CALCULATE BLADE SURFACE VELOCITIES I A c TER CONVERGENCE) 

C 

IF! TTER.NE.O) 63 TO 260 
00 250 K-l.KWX 

CALL SPLINE (SH(1,K),TT(1,K) ,MX,OELBTA,AC) 

A*XN 

no 240 1 * 1 , nx 

240 ABIf )*(RIT,K)*W*WAtI,K)*SRE7A(l,K))*f6.283186*R(I,K)/ A-TTI I»K> ) 
CALL SPLINE I SNI 1 ,K) ,43 ,HX,3R0M,AC) 

IF ISF4CT.LF. 1.0) GO TO 245 
A * SFACT*XN 
00 244 I=1,NX 

244 ABI I)*IR( I,<)*W«-WA{I,K1*SBETA{I,K))*t6.2B318S«R(I,K)/ A-TTI I,K)) 
CALL SPLINE I SNI 1,K) ,4S,NX, AO ,AC» 

245 00 250 I=1,MX 

BET AO = BETA C I ,X )-DEL3T4 ( I )/2, 

RET AT * BE1 A0+0ELBT4 IT) 

COSBO = COSIBETAO) 

C0S8T = CO S(BETAT) 

IFIZI I, K ) .LT .ZSPL IT) OROMII) * AO(I) 

WTRI I,K ) =Cf’SBO*C3SBT/(COSBO*COSBT)*I2.*HAII ,K) /CI)$8D*R| I,K)*W* 

II 8ETA0-BETAT ) /CBE'» A 1 1 , K >** 2+DRO MI I ) ) 

250 CONTINUE 
C 

0 END OF BLADE SURFACE VELOCITY CALCULATIONS 

C START CALCULATION 3F WEIGHT FLOW VS. 31 STANCE FROM HUB 
C 

260 DO 370 1*1, NX 
1ND*1 

00 270 K*1,KHX 
270 AC(K)*DN( I,K) 

GO TO 290 

280 WAI I, 1)*.5*WAI 1,1) 

290 00 300 K*2,XNX 
J*K-1 

HR*RI I,K )-R( I,J) 

HZ *7 1 1,X )-2l I, J) 

WAS*WA( I » J i*( l.*SA{ ', J)*HR+SCII ,J)*HZ)+SBt I »«l)*HR*SOf I »J)*HZ 
WASS*WA( I, J)+WAS*( S4 1 1 ,K)*HR*SC II ,F.) *HZ )*S8II ,K) *HR*SDI I ,K)*H" 

300 WA! I,K)*:wAS*WASS)/2. 

310 00 340 K*1,KNX 

TIP* l.-IWAI I,K)**2+2.*W*4LM-(W*R(1 ,K) )**2)/>./CP/TEHP 
IFIT1P.LT. .0) GO T3 280 

TPPIP* 1.- l2.*W*ALM-IW*RtI ,K) )**2)/2.7CP/TEHP 

" DENSTy*TlP**EAP3N*RH3-| T1P/TPP1 P) ♦•EXP0N*PL0SS/AR/TPP1P7TENP 
1 432. 17* SNI I,K)/SN(MXBL,K) 

PRM I, Kt^OENSTYPARPTlPPTEMP/BZ. 17/144. 

IFIZSI I ) .LE.7HI I ) ) SO TO 320 

PSI-ATANIIRSI D-RHII ))/;ZSII)-ZHII)) 1-1.5708 
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GO T3 330 

320 PSI=ATAN((ZH( I )-Z S( I ) )/(RS( II-RHU ))) 

330 WTHR'J=WA{ I , < )*C3ETA ( I ,K ) *C0S( PSI -AL ( I ,K) I 
A=XN 

!F(Z( I,K) .LT.Z SPLIT! A=SFACT*:CN 
C * 6.283186*R(I,K)-A«‘TT(l,K) 

340 AD(K)=DENSTY*WTHRU*C 

CALL INTGRLI AC( 1 ),A0( l) ,KMX,WTFL(1) ) 

IF ( ABS(WT-WTFL(KM J ) .LE . WTOLER) GO TO 350 
CALL CONTIN ( WA( I , 1 ) , WTFL ( KMX) , IND, I , WT) 

IF ( IND.NE .6 ) GO TO 290 
350 CALL SPLINT ( WTFL , AC , KMX ,BA ,KMX,AB) 

00 360 K =1 , KM X 

DELTA=ABS( AB ( K J-DNI I,KU 

ONI I,K ) = ( l.-CORFAC)*DM(I,K)«-CORFAC*AB(K) 

360 I FI DEL TA .GT. ERROR ) ERR OR =DELTA 
370 CONTINUE 

END OF LOOP - WEIGHT FLOW CALCULATION 

CALCULATE STREAMLINE COORDINATES FOR NEXT ITERATION 

00 380 K=2,KMXM1 
00 380 1=1, MX 

Z ( I, KI=DN( I,K I/ONI I»KNXI*(7$(I l-ZH( I ) l+ZHII ) 

3«0 R( I,KJ=DN( I, K ) /ON (I»KMX)*(RS(I J-RMI I))+RH(II 

IFC I ERRJR .GF.ERR0R1 ) . OR . I ERROR. LE. TOLER) I I TER* ITER-1 
IF( ITER.GT.O) GO TO 410 
WRITE (6,1100) 

00 400 K=l, KMX, NPRT 
WRITE (6,1080) K 
DO 390 1=1, MX 

AB( I )=( Z( I,K l-R ( I ,K) I /ROOT 
390 AC ( I ) =( Z ( I »K ) <-3 ( I , K » ) /ROOT 

CALL SPLINE ( AB , AC ,MX, AO ,CURV( 1 ,K) ) 

00 400 1=1, MX 

CURVI I, K )=CURV! I,K)/( l.+AD ( I ) **2)**1.5 

A=ON( 1,0*12. 

B= Z( I,K)*12. 

0= R( 1,0*12. 

400 WRITE (6,1110) A , B ,0, WA ( I ,K ) ,PRS (I , K) , WTR ( I ,K) , CURVI I ,K ) 
WRITE (6,1130) 

410 A=FRR0R*12. 

WRITE (6,1120) I TNO, A 
ITNO= ITN3 + 1 

IF (ITER.GE.O) GO TO 150 
WRITE (6,1140) 

K = (KMX+l)/2 
DO 440 1 = 1, M.1 

IF(ZS( I ) .LE. ZH( I ) ) GO TO 420 

PSl = ATAN((RS(I)-3H(I))/(ZS(I)-ZH(I)))-U5708 
GO TO 430 

420 PSI = ATAN((ZH(I)-ZS(I))/(RS(I)-RH(I))) 

430 AB( I ) = ( ON ( 1,2 )-OV(I,l >) *COS(PSI-AL( I ,2) I 
AC( I ) » (ONI I,K + l)-3N(I,K-l)*/2.*C0S(PSI-AL'.(fK)l 
A0( I ) = (ONI I,KMX)-DN(I, KMXM1) ) *COS(PSI-AL(I , KMX-1 1) 

A = A8( I )* 12 • 
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B = ACI I)*12. 

n = add )*12. 

444 WRITE (6, 1110) A,B,D 

CALL SPLINT I XT, THT4 , MTHTA ,Z 1 1 ,1 ) , MX,TNH ) 

CALL SPLINT 1 XT,THTA ,MTHTA ,Z (1 ,2) ,MX, THH) 

CALL SPLINT l XT ,THTA » MTHTA • Z( 1 »K ) ,MX, THM) 

CALL SPLINT IXT,THTA*MTHTA ,Z(l ,KNX ),MX,THS) 

PION * 3, 141 5927/XN 
TP ION * P10N*SFACT 
WRITE 16, 1150) 

DO 450 1=1, NX 

THHim = THHI I)-TPION>TTII,l)/?./RII ,1) 

THM1ID = THf*!D-TPlON*TTII,K)/2./RtI,K» 

THSH I) = THSID-TPION+TTII ,KMX )/2./RII,KMX ) 

THHKHID = THHID- PF ON-TT ( l , 1 ) /2. /R! I , l I 
THMXHI I ) = THM I II- PI ON-TT (I,K)/2./R(I«K) 

THSXHIII = THSID- PION-TT (I ,KNX )/2./RII,KMX ) 

THHXP(I) = THHID- PI ON*TT! I , l ) /2. /R (I ,1) 

THMKPII) = THHID- PION+TTII ,K)/2./RII ,K) 

THSXPII) = THSID- PICN+TTII ,KMX ) (?.• /Rll «KHX) 

THHXMXI I i= THHtD-TTlI.l )/2./RII,l ) 

THMXMXt 11= THHID-TTII,K )/2./RII,K ) 

THSKMXI II* THSI D-TTII , KMX ) /2. /RI I ,KMX) 

450 WRITE! 6, 11601 THH1CD .THHKHID ,THi)KP(D ,THHKHX( D,THM1 II) ,THHXHI I ) 
1,THMKP* D,THHKMX!D,THSim ,THSKHII ) ,T4SKP(I> ,THSKHX! I ) 

00 460 J = l»3 
1*1 


460 


1010 

1020 

1030 

1040 

1050 

1060 

1070 


K = 1 

IFIJ.EQ.2) K=lKNX*l)/2 
IFIJ.FQ.3) K=KMX 

TIP* l.-IWAI I,K)«-*2«2.*W*ALM-IW*RII ,KI )**2>/Z ./CP/TEMP 
DENSTY * TlP**EXPON*RHO 
C * 6.283186*RII,KI-XN*TTII,K)*SFACT 
WIDTH * rtB 

IFIJ.EQ.2) WIOTK * AC 
IFIJ.EQ.3) WIDTH * AD 
WH * 8AI 2 1 /OENSTY/C /WIDTH 
WTHFTA = ALM/R 1 I »K l-W*R 1 1 ,K) 

BETAIIJ) * ATANIWTHETA/WM) 

AAIJ) * 8 f IA II J 1*57. 29 577 
CONTINUE 

WRITE 16,1170) AA 

IFIBCDP.E0.2) CALL BCDUMPI SRW,BETAI (3) ) 

IFIRCDP.NE.l ) RETURN 

CALL BCDUMP I DNI 1 , 1 ) ,0N( 21 ,21 I ) 

CALL BCDUMP ( WAI 1 , 1 ) , WA ( 21 ,2D) 

CALL RCDUMP I 211,1), Z121,21>) 

CALL BCDUMP I Rll, 1), R(21,?1D 
RETURN 

FORMAT I4I5.6F10.4) 

FORMAT I8H1RUN NO . I 3 , 10X,25HI NPUT DATA CARO LISTING I 
FORMAT I7F10.4) 

FORMAT U0X24HBCD CAROS FOR DN,WA,2,R ) 

FORMAT I36HK STAG. SPFFO OF SOUND AT INLET * ,F9.2) 
FORMAT ( ///5X13HITERATI0N N0.I3) 

FORMAT I 1H 6X5HALPHA9X5HRC 9'.'4hSH 9X5HBETA 9X5HTT 


9X5H- 
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1 X5HSB 9X5HSr 9XSHS0 ) 

1080 FORMAT I 2X10HSTRFAMLINEI3) 

1090 FORMAT (9F14.61 

1100 FORMAT I1HL9X5H0N 15X5H7 15X53R 15X5HWA 1 5X5HPRESS 14X 3HM 

1TR14X3HRC J 

1110 \ URMAT 16F19.6.F18.6) 

1120 F3RMAT 1 18H t TERATION NO. I 3 , l OX, 2 '+HMAX. STREAMLINE CHANGE = , 

IF10.6 I 

1130 FORMAT ( LHJ ) 

1140 FORMAT I34H1 STREAMLINE SPACING ALONG N0RMAL//1 0 X5 HHIIB 14X5HMEAN 
1 14X6HSHR0UD1 

1150 FORMAT ? 1H1, 56X, i 7H81 ADE CUOROt NA TE S/1 9X.6HHU8 ,38X,6HNEAN , 

1 38X,6HSHR3UO/3I5X,IH1,10X,4HKHMX, 7X ,4-IKHPl , 7 X f 3HKMX , 3X ) » 

1160 FORMAT ( 12F11.4) 

1170 FORMAT I ///1KL ,10X.20HINLE1 ANGLES - HU8tF7.2,8H, MEANF7.2, 10H 

1, SHR3U0F7 . 2 ) 

ENO 
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AIBFTC SPLBLO DFCK 

SUBROUTINE SPLBtO(HU3 ,HEAN,SH e .OUO,MOB) 

COMMON SRW, <MXF, MX,WT,XN , SAM, /iR, TYPE, BCOP, TEMP, AIM, RHO, TOLER, 

1 PLOSS.NP1T, ITERl,3ETIN,WTOI.ER,TMHl ,T-HKH,THHKP*THHKMX»THM1 , 

2 THKKH»THMKP«THMKMX» THS1 »THSKH»TMSKP»THSKMX»Zl ,R1 , AB, AO, AE, 

3 RIIMO, MXBL »P f ON, W, 3ETAI ,DTOM,CURV 

DIMENSION Z(21),R(21),ON(21),SM(21) ,BAI22> ,AB (22 ) , ACI22 », AL C 21) , 
lRC(2l),r.AL(21),SAL(21),PRS(2 2 ,WTFL (22 ) ,WTHETA ( 21 ,22) ,DWOM( 21 , 22 ) * 
2THETA( 21,22) ,WA( 21,22) ,BFTA(2*,22) ,S8ETA(21 ,22) ,CBETA(21 ,22 ) , 

3$A( 21, 22 ), SB (21, 22), C UR V( 21, 22) ,OT")M(2l ,22) 

DIMENSION Zl(21,21),Rl(2l,21t , THHi (21) «THHKH(21 ) »THHKP( 21), 

1 THHKMXC 21 ) , THN1 ( 21 ) »THMKH(21) ,HMKP(2l I ,TMMKMX(21 ) ,THS1 (21 1, 

2 THS< H( 21 ) »THSKP ( 21 ) , THSKMX(21 ) 

DIMENSION THAI (23,4) , SMAL (23) ,AO (22 ) ,AE(?2 ) ,DENSTV( 22 ) , BETAK 3) 

INTEGER RUNO, TYPE ,BCTP ,SRM ,RR ,HUB ,SHR(JUD 

LOGICAL PRINT 

IFIMOB.EQ.2) GO TO 10 

RUN0=0 

IF(MOB.EQ.A) GO TO 10 

S REAO (5, 1010)MX,KNX»MX$P,RR ,W,WT ,XN,GAM,AR 

ITNO * 1 

RUNO*RUNO+i 

WRITF (6,1020) RUNO 

WRITE (6, 1010 )MX, KMX, MXSP,RR ,W,MT ,XN,GAM,AR 

REAO (5, 1010 ) TYPE »3S3P »SRW » MXBL, TEMP f ALN, RHO, TOLER, PL OSS » Oil DM1 

MRITE(6, 1010)TYPE,B:DP, SRW, MXBL, TEMP, ALM,RNO, TOLER, PLOSS, ON DM1 
PL0S$*PL0SS*144. 

READ (3,1010) NULL, NPRT, ITER, NULL, 9ETIN,HT0LER,C0RrAr. 

WRITE (6, 1010) NULL ,N?RT,ITER ,NULL,BETI N,MTOLER,CORFAC 

BETIN a BETIN/57. 29577 

KHMX - KMX/2 

KHP 1 * KHMX+1 

KHP2 * <HMX*2 

REAO (5,1030) (THETA (1,1), 1*1, MX) 

WRITS(6, 1030) (THETA (1,1), 1*1, MX) 

READ ( 5, 1030)1 THETA 1 1 ,KMX) ,1*1 , MX) 

WRITE! 6, 1030) (THETA (I, KMX) ,1*1, MX) 

READ (5,1030) ( THETA( I ,KHMX) ,1 *1 ,MXSP) 

WRITE! 6, 1030) ( THETA ( I ,KHMX) ,1 *1 ,MXSP) 

READ (5,1030) ( THETA ( I ,KHP1) , 1*1 ,MXSP) 

WITITE(6, 1030 ) ( THE T A ( I , KHP 1) ,I=1,MXSP) 

IF(RR.EQ.l) GO TO 90 
REAO (5, 1030)12(1), 1*1, MX) 

WRITE! 6, 1030) (2 ( I ),!*1,MX) 

REAO ( 5, 1030)(R( I )»I=1»MX) 

WRITE(6,1030)(R( I), 1=1, MX) 

REAO ( 5, 1030) I ON ( I), 1*1, MX) 

WRITE! 6, 1030 ) (0N( I), 1*1, MX) 

DO 9 1*1, MX 

ZC i)=z; U212. 

R( I )*R( I ) / 2. 

9 ONI I ) =0N( 1 1/12* 

GO T3 21 
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10 IF( MO B.EQ » 4) CALL 3CREAD( SRW,8ETAI(3) ) 

MX = MX+1 

MOB = 2 

WT = WT/FLOATIKMXF-l ) 

WTOLER = WTOLER/FLOAT(KMXF-1 ) 

MXBL = MXBL+1 

11 READ (5, 1010)KMX»MXSP,RR t TYPE , CORF AC .TOLER, THI KMX 
RUNO * RUNO+l 

IF C <HUB.EQ.OJ .AND.(MEAN»EQ.O) . AND. ( SHROUD. EQ. 0) ) RETURN 
© WRITE (6,1020) RUNO 

•g} WRITE(6, 1010)KMX,MXSP,RR , TYPE .CORFAC .TOLER, THI KMX 

to ITN3 = 1 

W ITER = I TER l 

KHMX = KMX/2 
K i IP 1 = KHMX+1 
KHP2 * KHMX+2 

READ (5,1030) Z( 1 ) ,R ( 1) ,DN ( 1 ) 

WRITE! 6, 1030) Z( 1),R( 1) ,0N(1) 

zm * zm/i 2 . 

R( 1 ) = R ( 11/12* 

DNt 1) = DN(1)/12. 

DWDM l = W+ALM/RC l)**2 
|F( HUB.EO.O) GO TO 14 
HUB * 0 
DO 12 (=2, MX 
J = 1-1 

THETA! 1, 1) = THHKJ) 

THETA( I, KMX) = THHKMX! J ) 

Z( I) « Zl( J, 1) 

RID * R It J, 1 ) 

12 DN( I) = AB(J ) 

DO 13 !=2»MXSP 
J * 1-1 

THETA ( I.KHMX) * THHKH(J) 

13 THFTA( I,<HP1) * THHKP(J) 

BET IN * BETAK l) 

GO TO 20 

14 I F ( MFAN.EQ .0 ) GO TO 17 
MEAN » 0 

K * ( KMXF + D/2 
DO 15 1=2, MX 
J = 1-1 

THET A ( I, 1) = THMl(J) 

THETA ( I, KMX) = THMKMX! ") 

Z( I) = 2 1 ( J , K ) 

R( I) = R 1 ( J , K ) 

15 DN( I ) = AD( J ) 

DO 16 I=2»MXSP 
J = 1-1 

THETA( I.KHMX ) = THMKH(J) 

16 THETAC I, KHP1 ) = THMKP(J) 

BETH = BETA 1(2) 

GO TO 20 

17 IF ( SHROUD. EQ.O) RETURN 
SHROUD = 0 
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00 10 1=2, MX 
J = 1-1 

THETAU, 1) = THS1CJ1 
THFTAll.KMX) = THSKMX(J) 

Z(I) = Z 11 J.KMXF ) 

R(I) = Rll J.KMXF) 

18 0N( I) = AEIJ ) 

00 IP 1*2, MXSP 
J * 1-1 

THFTAU.KHMX) = THSKHl J) 

19 THETAl I.KHPl ) = THSKPIJ) 

BETIN = BETA 1(31 

20 TANBTO = SINl BETIN)/COSlBFTlN) 

TANBT = TAN8T0*( ALM-W*Rl 1 )**2) / l ALM— W*R( 2 1 **2 ) *DN( 1 )/ON(2 ) 

DTDM1 = TANBT/RU) 

OELTHT = TANRTO/(W*Rl **2-ALM)* ( ALM*ALOG (R(l)/R(2) )*W/2.*(R( 2)**2 
1 -R( 1 )**2) '* ( ONI 1)*0N( 2) ) /2» /ON (2 ) 

THETAl 1,KMX) * THIKMX*OELTHT 
THET A ( 1, <HP 1 1 = THETAl l,KMX)-Pl ON 
THETAl 1, KHMX) = THETAl 1 ,KHP1 ) 

THETAll, 1) = THETAl 1,KHMX)-PI ON 
NULL = 0 

WRITE 16, 1 010 )MX»KMX »MXSP»RR ,W , WT ,XN,GAM,AR 
PL0SS1 = PL3SS/144. 

WRITE 16, 1 0 10 )TYPE »9«DP»SRW»MXBL,TEMP s ALM ,RHO» TOLER, PLOSSl.DWDMl 
BETIN 1 = BET IN* 5 7. 295 77 

WRITE16, 1010) NULL,NPRT , I TER .NULL ,BETI N1 ,WTOLER»CORFAC 
WRITE16, 1030) l THETA 1 1,1), 1=1, MX) 

WRITEl 6, 1030) (THETA 1 1, KMX) ,1*1 , MX) 

WRITE! 6, 1030) (THETAl I, KHMX) ,1*1, MXSP) 

WRITFl 6, 1030) ( THETAl I, KHP1) ,1*1 ,MXSP) 

WRITEl 6, 1030) IZ (I), 1*1, MX) 

WRITEl 6, 1030) IRl I ),I*1,MX) 

WPITE16, 1030) (DNl I), 1*1, MX) 

21 IFl TTPE.EQ . 1 ) GO TO 30 

WA( 1,1 )=WT/RHO /ONI l)/Rll)/XN/( THETAll, KMX) -THLT All, 1) ) 

DO 23 1*1, MX 

IFl ( I .EQ. D.OR.l l.GT.MXSP) ) THFTA (I ,KHMXJ * ( THETAl I ,KMX) +THET A 

1 ( I, 1) )/2. 

IFl ( I. EO. D.OR.l l.GT.MXSP)) THETAll ,KHP15 * THETAl I , KHMX) 

00 22 K*l,KHMX 

THETAl !,<)'» FL0AT(K-1)/FL0ATIKHMX-1)* (THETA II , KHMX) -THET Al 1,1)) 

1 «• THETAl 1,1) 

22 WA( I, K ) =WA 11,1) 

00 23 K*KHPl ,KMX 

THETAl I , K ) * FLOA T( K-KHPi ) /F LOATl KMX-KHP1 ) *( THET A( I ,KMX )- 
1 THETAl I, KHP1) )+THETAU ,KHP1) 

23 WAl I , K > * WAll.l) 

REAO l 5, 1010) NFXT 


24 


WRITE(6»1010) NEXT 
IFINEXT.EO.O) GO TO 25 

REAO (5, 1021) li.Kl.THFTAUl ,K1) » 1 2, K? »THtTA( I 2 ,K2) ,I3,K3, 
1 K3)»I4,K4,THETA(I4»K4) . I B»K5 ,THE TA II 5 , K5) .NEXT 


THETAl 13, 


WRITEl 6, 102 1)11, Kl, THETA (II ,K1) ,12 ,K2 ,THFT Al 1 2 ,K?) , 13 ,K3, THETAl 13, 
1 K3),T4,K4, THETAl 14, K4) ,?5,K5,THETA(I5,K5) ,NEXT 
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GO TO 24 
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ENO OF INPUT STATEMENTS 

25 WRITE (6, 102?) 

00 26 K=1,KWX 

26 WRITF (6, 1030) I THETA ( I , K) ,1 =1 , MX) 

GO TO 45 

30 CALL BCREAO ( THETA ( 1, 1 1 ,THETA(2i ,21 ) ) 

CALL BCREAD ( WA (1,1), WA(21,71)I 

WRITE (6,1040) 

45 CONTINUE 

CP=AR*GAM/(GAM-l. ) 

Cl = SQRT(GAM*AR*TEMP) 

WRITE (6,1050) Cl 
KMXMl * KMX-1 
MXP2 * MX + 2 
MXSP1 * MXSP+1 
WTH = WT/2. 

EXPON * 1 • /( GAM- 1 . ) 
DT0M1=SIN(BETIN)/C0S(BFTIN)/R{1) 

IFOOB.EQ.2) 0TDM1 = TANBT/R ( 1 ) 

ROOT=SQRT( 2. ) 

SMI 1 1 * -SQR T( (Z(2)-Z(l>)**2+(R(2)-R (111**2) 

00 60 K = l, KMXMl 

60 BA( K ) =FLOAT( K- 1 ) *WT/FL0AT( KMXM1-1 ) 

CALCULATE ALPHA AND SM 

DO 70 1=1, MX 
AB(I)=(Z(I)-R(I))/ROOT 
70 AC(I»*(Z( I )+R 1 1 ) ) /ROOT 

CALL SPL INE( AB,AC,MX,AL,RC) 

DO 80 1=1, MX 

AL ( I ) =ATAN ( AL ( I ) >-.785398 
CAL ( I ) =C3S( AL ( I ) ) 

80 SAL ( I ) =S lNj( AL ( I ) ) 

DO 85 I=2»MX 
J=I-l 

85 SM( I)=SM( JUSQRT( ( Z ( I )-Z( J) ) **2 + ( R( I) -R( J) ) **?) 
SMAL13) = SW( 2>+.l*(SM(3)-SM(2> ) 

SMAL( 4) = SM(2)+.5*(SM(3)-SM(21 ) 

90 ERROR = 1000. 

CALCULATE BETA ON BLADE SURFACES 
BEGINNING OF LOOP FOR ITERATIONS 

91 DO 97 J = 1 , 4 
K = 1 

IF( J.'3.2) K=KHMX 

I F( J . EQ . 3 ) K=KHP1 

IF( J.EQ.4) <=KMX 

00 93 1 = 1,2 

THAL( I, J > = THETA ( I ,K ) 

93 SMAL( I) = SM ( I ) 
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DO 94 r = 3, 4 

94 THAI ( If J ) = THETA ( 2 »K ) + C SMALC I )— SMC 2 ) ) /C SMC3 ) -SMC 2) )*CTHET At 3, X t 
l -THETAt 2»K ) ) 

DO 95 I=5»MXP2 
IM2 = 1-2 

THALC I,J ) = THETAC IM2,K> 

95 SM AL ( I) * SNt IM2> 

CALL SPL INt 2t SMAL , THAL < 1 J ) ,DTDM1 ,M:tP2 ,DTDMt l ,K) ,AB) 

DO 96 1=1,2 

CURVtI.O = ABCt )/12./C1.+CDTDMtI,K)/12.l**2k**1.5 
BETA! I,K ) = ATANCRtt)*OTOMCI,Kn 
SBETAC I, < ) = SINCBETA 1 1 ,K) ) 

C8FTA( I,K ) = CO St BFTA t l ,K) ) 

96 CONTINUE 

DO 97 1=3, MX 
IP2 = 1+2 

CURVCI.M = ABCIP2I '12./tl.+IDTDMCIP2,K)/12.)**2)**1.5 
BETA! 1,0 = ATANCRt I)*DT0M(IP2,K) > 

SB ET All, XI = SINCBETA ( I ,KI I 
CBETAII.O = COS( BFTA ( I ,K I ) 

97 CONTCNUt 

PRINT = ( ITER.LE.OI.OR.CITNO.LE.NCI 
IF I PRINT) HR ITE J 6, 10601 1 TNO 
ERROR l=E«ROR 
ERROR =0. 

START CALCULATION OF PARAMETERS 

LAST = KHMX-1 
DO 100 K=2*L AST 

CALL SPL IN2( SM, THETA Cl,K)»DTDMl,MX,DTDMtl*lO»A8l 
DO 100 1=1, NX 

CUR VC I ,< ) = AB<I)/l2./tl.+tOTDMCI,K)/l2.)**2I**1.5 
BETAC I,O=ATANCRCI>*DT0MCI,K)> 

SBETAC I,K )=SINCBFTACI ,K) ) 

- CBETAC I,K)=COSCBETAU,KJI 

100 CONTINUE 

DO 101 K=K HP2,KMXM1 

CALL SPL IN2CSM, THETA C1,KI,0TDM1, MX, DTOMl 1 , K» , AB) 

DO 101 1=1, NX 

CURV(I,<> = ABCI 1/12./C l. + COTDM CI,K)/12.) **2 ) **1 .,5 
BETAC I,X)=ATANCRC n*DTDMCI,KH 
SBETAC I,K)=SINCBETACI,Kn 
CBETAC !,K)*CO$CBETAtI,Kn 

101 CONTINUE 

DO 110 K*1,KMX 
DO 105 1=1, MX 

MTHETACI,K»=MACI,K)*SBFTACT,K) 

105 CONTINUE 

CALL SPL IN?( SM,WTHE.TA ( 1,KI *DWDM1 ,MX ,0H0M<1 ,K» , AC» 

DO 110 l«l,MX 

SAC I, K I ■ SAI.C I)*SBETACI,K»*CBFTACI ,KI 

SBC I • K » = CBETAC I,K)*RCn*C2.*W*SALCn+C/H0MCI,tU) 

110 CONTINUE 

ENO OF PARAMETER CALCULATION 
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C CALCULATE VELOCITY 0 I STR 13 UTI ON t CHECK CONTINUITY 
C 

00 ?00 1=1, MX 
TN0=1 
IND2 = 1 
GO TO 130 

120 WAI I, l)=.5*WAl 1,11 
GO TO 130 

125 WAI I , 1 ) =2.*WA 11,1) 

GO TO 130 

126 WAll.KHPi) * .5*WA II , KHPl) 

GO TO 130 

127 WAI I » K HP l ) = 2.*WA(I ,KHP1) 

130 IFHI.GT.1).AN0.1I.LF.MXSP)) 50 TO 142 
00 140 K=2»KMX 
J * <-l 

IFU.EQ.KHPl) WAII.K) = WAII,J) 

IFU.EQ.KHPl) GO TO 140 
HT=THETAl l,K)-THETA!l ,J) 

WAS = HA I I» J )+( WAI l,J)*SAII,J)«-S8II » J) )*HT 
WASS = WAI l,J)*IWAS*SAII,K)«-SBU,K))*HT 
WAU,K) = I WAS+WA SS) /2. 

140 CONTINUE 
GO TO 148 

142 00 144 K=2,KHMX 
J=K-1 

HTsTHETAt I ,K )-THETA I I , J) 

WAS = WAI I,J)+IWAII,J)*SAII,J)*SBII,J))*WT 
WASS * WAI I,J)+lWA$*SAII,K)+SBII,K) !*HT 
WAI 1 »K ) * IWAS+WASS)/?. 

144 CONTINUE 

90 146 K=KHP2»KWX 
J = K-l 

HT*THETAI I ,K )-THETA l I , J > 

WAS « WA(I,J)*IWAII»J)*SAII»J)*SB(I»J))*HT 
WASS « WAI I,J)+IWAS*SAU,K)+SBII,K) )*HT 
WAI It K I = I WAS+WA SS) /?. 

146 CONTINUE 
148 CONTINUE 

DO 150 K=1,KMX 

TIP* 1.-IWAI 1,K)**2+2.*W*ALM-IW*RII) 1**2) /2 ./CP/TFM* 

1FI I I 1P.LT..0) .AND.II.GT.D.AND. II.LE.MXSP). A NO. (K.GE.KHPl ) ) 
l GO TO 126 
IFITIP.LT..O) GO TO 120 

TPt>lP= 1.- 1 2.*W*ALM-IW*R(l I )**2) fZ. /CP/TEMP 

DENSTYIK) * T1P**EXP0N*RH0-I Tl P/TPP1P) **EXPON*PLOSS'AR/TPPlP/TEHP 
l *32. 17*SMI I ) /SMINX8L) 

PRSt < ) = OENSTYI K )*AR* TIP* TEMP/3 2. 17 /144. 

WM=WA I !,K)*CBETAI I,K) 

ABU) = DENSTYIK )*WM*9NI I)*RII l*XN 
150 ACU)«THETAII,K) 

CALL INTGRLIAC,AB,KHNX,WTFL) 

IFIWTFLtKHMX) .LE..O) GO TO 125 
CALL INTGRL I AC IKHP1 ) fABIKHPl) ,KHMX,WTFLIKHP1) ) 
IFIIWTFUKMX).LE..0). AND. II. GT.l). AND. II.LE.MXSP)) GO TO 127 
IFIWTFLI KMX).LF..O) GO TO 125 
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IF( ( I.GT.l).AND.(I.LE.MXSP)) GO TO 155 
DO 152 K=KHP ltKMXHl 
ACIK) = ACIK+l) 

152 WTFLI K ) = WTFL (KHMX) +WTFL I Kfl ) 

IF(ABS(WT-WTFL(KMXMl)).LE.WTOLFR) 30 TO 150 
CALL CONTIN(WA(I,l),WTFL(KMXMl) ,INO,I ,WT) 

IF! IN0.NE.6JG0 TO 130 
GO TO 160 

155 IF I ABSI W TH-WTFL (KHMX)).LE.WTOLFR) IN0*6 
IF ( A8S( WTH-HTFL (KMX ) > .LE.WTOLFR) IN02=6 
I F ( INO. NE. 6) CALL CONTI Nl MA (I , l ) ,WTFL(KHMX> , INO, I ,WTH) 

IF! IND2.NE.6f CALL CON TAL ( WA (1 »KHP1 ) » WTFLI KMX) »IND2 » I ,WTH ) 

IF(( IND.VF. 6 ). 0 R.{INO 2 .NE. 6 n GO TO 130 
160 IF(( I .GT.U.ANO.I I.LE.MXSPi) GO TO 165 
CALL SFlINT (WTFL ,AC ,KMXM1 ,9A,KMXMl ,AB) 

GO TO 166 

165 CALL SPL INT( WTFL »AC » KHMX,BA , KHMX, AB ) 

CALL SPLINT(HTFLIKHPl) ,AC(KHP1) , KHMX, BAI1) ,KHMX,A8 (KHPl)) 

166 CONTINUE 

IFf(I.GT.l).AND.II.LE.MXSP)) GO TO 168 
00 167 K=1 f KHMX 
KA =* KMX-K+1 
J = KA-1 

167 ABIKA) * AB(J) 

168 00 170 K = 1 , KMX 
DFLTA=ABS(AB(K I- THETA 1 1 5 K) ) 

170 IF! DEL TA.GT. ERROR) ERROR =DELTA 
IF( .NOT. PRINT) GO TO 178 
A=$M( I >*12. 

C*AL( I 1*57.29577 
0 * R ( I )*12. 

E * ZI 1)*12. 

F = ONI I )* 12. 

WRITE (6,10801 I,A,C,0,E,F 
WRITE ( 6, 1 070 ) 

00 175 K=1,KHMX»NPRT 
B*BETA( !,:< >*57.29577 
C * W THE TAJ I ,K >+W*R ( I > 

WM * MAI I ,K >*C5ETA 1 1 , K ) 

V = SORT I C**2+WM**2 ) 

OWOT - WA(I f K)*SMI,K) + SB(l»K) 

WRITE (6,1090) THETA II ,KI »CURV(I ,KI ,B,WA(I,K| , WTHETAI I,K),C,MM, Y, 
1 PRSIO,OENSTY(K>,DTOM(I,K) ,DWOM( I ,K) , S A 1 1 , K> ,SBI I ,K) ,DWDT 

175 CONTINUE 
WRITE C 6, 1095 > 

DO 176 K *K HP 1 » KM X » NPR T 
8° BETA! I ,K 1*57.29577 
C = WTHETAI l,KHW*R(I ) 

MM = MM l,K)*CBETA,I,KJ 

V - SORT! C**2<WM**2 ) 

OMOT « MAI I ,K )*SA( I ,K l+SBI I ,K) 

WRITE (6,1090) THETA ( I ,KI ,CURV(I,K) ,B ,MA(I ,K> , WTHETAI I ,K» ,C,WM, /, 
l PRSf X ) » OENSTYI K), DTD Mil, K) ,OWOM( t,K),SA(I,K),SB(I,K) ,DWOT 

176 CONTINUE 
WRITE (6,1095) 

178 DO IS'' K»2,KMXMl 


u o u 


180 THETAJ I* < I =( l.-CORFAC )*THFTA ( I ,K)f;0RFAC*AB(X) 

200 CONTINUE 

END OF VELOCITY CALCULATIONS. 

WRITE (6,1120) !TNO, ERROR 
IF( ITER.LE.O) GO TO 230 

IF( ( ERROR. GE. ERROR 1). OR. (ERR OR. LE. TOLER) ) ITER*ITER-1 
ITNO= ITN2 + 1 
GO TO 9i 

230 IF(BCOP.NE^l) GO TO 240 

CALL BCQUMP ( THETA ( 1, l ) ,THF (21 ,21)) 

.CALL BCDUMP ( WA (1,1), WA (21,21)) 

240 IF(M3B.EQ.2) GO TO 11 
GO TO 5 

1010 FORMAT (4I5,6FK.4) 

1020 FORMAT (8HIRUN NO. 13, 10X,25HINPUT DATA CARD LISTir ) 

1021 FORMAT < 5 ( 21 2,F8. 5 ) , 1 1 ) 

1023 FORMAT ( 32H THETA-CALCULATEO AND/OR INPUT ) 

1030 FORMAT (7F10.4) 

1040 FORMAT (10X2';HBC0 CAROS FOR THETA, WA ) 

1050 FORMAT ( 36HK STAG. SPEED OF SOUND AT INLET » ,F9.2,//( 

1060 FORMAT ( ///5X13HI TERA TI ON NO. 13) 

1070 FORMAT ( 132HK THFTA T-CURV SETA WA WTHET A VTHFTA 
1MH V PR S 9ENSTY DTDM DWDM SA SB 

20T ) 

1080 FORMAT ( 2X16HQUASI-0RTH0G0NALI 3,6X,4HSM »,F7.4,9H ALPHA *, 7.2 
1 5H R *,F7.4,5H l =,F7.4,6H ON *,F7.4) 

1090 FORMAT ( 1X,F9.4, 2FP.2,5F8. 1 ,F7,2,F9.5 ,F7.2,F8.0,F8>4,2F9. 1) 

1095 FORMAT C 1H ) 

1120 FORMAT ( 18HJ ITERATION NO. 1 3 ,10X,24HMAX. STREAMLINE CHANGE 
1F10.6) 

1200 FORMAT! //10X7HNORMAL 14) 

1210 F0RMAT(7F18.6) 

ENO 
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SI8FTC SPLINT 0FCK 

SU8R3UT INF SPLINT < X , Y,N, l , MAX , YI NT) 

DIMENSION X< 501 , Y( 50 ) , S( 50) ,A(50) ,3 15 0) ,C (50) , F (50 I ,W(50) ,SB150 ) , 
1GI 50>,EM( 50),Z< 50 ) , YI NT( 50) 

C0MM3N 0 
INTEGER Q 
00 10 1=2, N 
10 SUhXUHIM) 

NO=N~l 

00 ZU 1=2, N3 
AlI)=Sin/6.0 
BUh(S( I)+SIIU)l/3.0 

cc n=s( i + 1 )/6.o 

20 fcimiyi iHi-y(:))/$(Hu-(Y(n-y(i-in/su> 

A(N)=-.5 
&< 11*1.0 
B( N)= 1 *0 
CU)*-.5 
FC 11*0.0 
Ff NI=0.0 
HI 1 >*811) 

SBI lKm/Kll) 

Gf 1 )*0.0 
00 30 1*2, N 

Hm=8{fNU)*S9(HI 

$Bm*cm/H(n 

30 GII)*CFm-AI I)*GU-l)l/wm 
6M(N)=G<N) 

00 40 1*2, N 
K=N+l-I 

40 EM(K )=G( X l-SBCK )*EM( K+ll 

oo oo i=i, M ax 

X=2 

!FfZm-X(l>) 60,50,70 
50 YINTI I)=Y( 1) 

GO TO 90 

60 IF(ZCU.LT.( 1.1*X( I)-. t*X(2)MWRITE (6,1000>Z(I> 

GO T3 85 

1000 FORMAT ( 17H OUT OF RANGE Z =F10.6) 

65 IFCZ< D.GT.I 1. 1*X(N)~. l^Xf N-l ) )) WRITE (S,lD33IZ(n 
K*N 

GO TO 85 

70 IFIZI I)-X(K>) 85,75,80 
75 VINTC n*Y<K) 

GO TO 90 
80 K*K*1 

IFIK-N) 70v 70, 65 

85 Y INTI 1 1 * E4(K-nMXm-Z<in**3/6./$(K) + EM<K)*<Zm-X(K-l))**3/6. 
i/SliOmm/ <K)-FKCK)*$(K)/6.l*(Zm-X(K-l)H>mR-l)/S|K)-EMtK-il 
2*S(K)/6.>*fX(K)-ZU>) 

90 CONTINUE 

MXA * MAXOI N,MAX ) 

1FIQ.E0.16) WRITE (6, 1010) N,MA X, ( X( I) , Y( I) , Z ( I I , YINTI I),!*i,NXA) 

1010 FORMAT (2X21HN0. OF POINTS GIVEN *,!3,30Ht NO. OF INTERPOLATED POI 
INTS *, 13, /10X5HX 15X5HY 12X11 HX-INTERP0L.9XHHY~INTERP0L./f 4 


2E2C.8 } ) 
100 RETURN 
ENO 



E-3Z98 


21 


$IBFTC SPLDFR OECK 

SUBROUTINE SPLOERt X , Y,N,Z ,MAk,D YDX) 

DIMENSION XI 50) ,Yt SO), SI 50) ,AI50) ,8150) ,CI53) ,F 153 ) ,WI 50) , S8I 50 » , 

1G( 50); FMt 50) ,21 50) »OYDX( 50) 

00 10 1=2, N 
10 SI I I =X( l . -X( I- 1 ) 

N0=N-1 

00 20 1=2, NO 
A( I)*SI D/6.0 
B( I)=(S( I»+$I I + 1D/3.0 
C( I) = S( I "D/6.0 

20 Fm*im+l)-YU))/S(’ + l! lYm-Yll-lM/SII) 

A(N)=-.5 
BI 1) = 1.0 
B(N)=1.0 
Cl l)=- .5 
FI 11=0.0 
F(N)=0.0 
Wl l)=8tl) 

SB! 1)=CI 1 S/HI i) 

61 1 )=0.0 
00 30 1 = 2, N 

Wl I )=BI D-AI £ »*SB 1 1— 1 * 

SB! I ) =Cl 1 )/Wt I ) 

30 GI I > = t FI I )-A I I )*G 1 1- 1 ) ) /WI I ) 

EHIN)*GtN> 

00 40 1=2, N 
X=N*1-1 

40 EMIK) =GI< )-SB # K )*EMIK*11 
00 90 1=1, MAX 
K*2 

IF1ZI II-XIIH 60,70,70 
60 WRITE 16,1000)211) 

1000 FORMAT I 17H 0U1 OF BLADE 2 =F10.6) 

GO TO 85 " 

65 WRITE 16, 1000)ZI I ) 

K=N 

GO TO 85 

70 IF1ZI I)-X(K) ) 85,85,80 
80 K*K«-1 

IFIK-N) 70,70,65 ,Jf. 

85 DYDXI I)=-EMIK-1)*(XIK)-ZII ) )**2 /2. O/SI K) ♦EMIK) 41X1 K-l )-ZI D )**2/2. . r i 

10/St <)♦! Y( K J-YIK- 1) )/SlK)-lEMIK)-EHIK-l))4SIK)/6.0 V»- 

90 CONTINUE 
100 RETURN 
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SIBFTC CONT IN DECK 

SUBROUTINE CONTI N ( WA ,WTFL ,IND, ! ,WT) 

DIMENSION SPEED! 3), WEIGHT! 3) 

135 GO T3 (140,150,210,270,370) , IND 
140 SPEED! 1) = WA 
WEIGHT! 1) = WTFL 
DEI T A * WT/WTFL*WA-WA 

IF( ABS! DELTA ) .GT • 100. ) DELTA = SIGNUOO. , DELTA) 
WA = OELTA+WA 
IND * 2 
RETURN 

150 IF ( i WTFL-WEIGHT(l) )/( WA-SPEEOI I) ) ) 180,190,160 
160 SPEED! 21 = WA 

DELTA = (WT-WTFL)/(WTFL-WE1GHT(1) )*(WA-SPEED(l) ) 
IF(ABS(0£LTA).GT.100.) DELTA = SIGNUOO. , DELTA) 
WA = DEL TA ♦WA 
166 SPEED! 1) = SPEED! 2) 

WFIGHT(l) = WTFL 
RETURN 

170 WRITF (6,1000) I, WTFL 
IND = 6 
RETURN 
180 IND = 3 

IF (WTFL.GE.WT) GO TO 140 
IF (SPEED! D-WA) 190,200,200 
190 SPEED! 2) = SPEED! 1) 

SPEED! 1) = 2 . 0* SP FED ( 1) - WA 
SPEED! 3) --- WA 
WEIGHT! 2) = WEIGHT(l) 

WEIGHT! 3) * WTFL 
WA = SPEED! 1) 

RETURN 

200 SPEED! 2) = WA 

SPEED! 3) = SPEED! 1) 

SPEED!)) = 2.0*WA-SPEE0( l) 

WEIGHT! 2) = WTFL 
WEIGHT! 3) = WEIGHT! 1 ) 

WA * SPEED! 1) 

RETURN 

210 WEIGHT! 1) = WTFL 

IF (WTFL.GE.WT) GO TO 140 
IF ( WEIGHT ( 1 )-WE IGHT! 2) ) 230,380,220 
2 20 WEIGHT! 3) = WEIGHT(2) 

WEIGHT! 2) = WEIGHT! 1) 

SPEED! 3) = SPEED! 7) 

SPEED! 2) * SPEED ( l ) 

SPEED! 1) = 2.0*SPEFD(2)-SPEED(3) 

WA * SPEED! 1) 

RETURN 

230 IF (SPEED! 3)- SPEED! lJ-10.0) 170,170,240 
240 IND « 4 

245 IF (WEI GHT ( 3 ) - WE I GHT ( l ) ) 260,260,250 
250 WA = ( SPEED! I ) ♦SPEED! 2) ) /2. 0 
RETURN 

260 WA » (SPEEOm+SPEFO! 2) )/2.0 
RETURN 
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270 IF 1SPEED! 3)-SPEED!l*-l0.0) 170,170,280 
280 IP ( 1,'TFL- WEIGHT! 21 I 320,350,290 
290 IF ( W A-SPEED! 2 ) ) 310,300,300 
300 SPEED! 1) = SPEED ( 2 * 

SPEEDC2) = WA 
WFIGHTU* = WEIGHT! 2) 

WEIGHT! 2) = WIFI 
GO T3 24 5 

310 SPEED! 31 = SPEED! 2) 

SPEED! 2) = WA 
WEIGHT! 3) = WEIGHT ! 21 
WEIGHT! 2) = WTFL 
GO TO 245 

320 IF I WA-SPEFD! 2 ) ) 340,330,330 
330 WEIGHT! 3) = WTFL 
SPEED! 3) = WA 
GO TO 245 

340 WEIGHT! 1 1 = WTFL 
SPEED! 1) = WA 
GO TO 245 
350 INO = 5 

IF ! WA-SPEEO! 2 1 1 380,360,360 
360 SPEED!!* * SPEED! 21 
WEIGHT! II = WEIGHT! 2* 

SPEE0I2) = I SPEED! 1 l+SPESD ! 31 1 /2« 0 
WA = SPEED! 2 1 
RETURN 
270 INO * 4 

WEIGHT! 2) = WTFL 
WA =* ( SP EEDi 1 ) +SPEED! 2) * /2. 0 
RETURN 
380 IND = 5 

390 WEIGHT! 3) = WEIGHT! 21 
SPEED! 3) = SPEED! 2* 

SP FED! 2) = ! SPEED! 1 X-SPEEDI 3*1/2. 

WA = SPEED! 2 * 

RETURN 

1000 FORMAT I/12H FIXED LINE I2,12H, MAX WT * FID. 61 
END 
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MBFTC SPIIN2 DECK 

SUBROUTINE SPUN 2 ( X, Y , YIP , N , SLOPE ,E M) 

DIMENSION X{ 50) , Y ( 50 '» » SI 50) ,A(50) f B(50> ,C (53 ) ,F ( 5 j } »«f 50 ) ,S BI 50 ) , 
1G< 50 ) » EMI 50) ,SL0PE(50) 

COMMON Q 
INTEGER 0 
DO 10 1=2, N 
10 St I)=X( I )-X( I-i) 

N0=N-1 

00 20 1=2, N3 
A( I ) = SI l)/6. 

B( I) = ( S( I1+MH1) )/3. 

Cl I) = S( I + D/6. 

20 FI I) = (Y( lUI-YIID/SIIUI-IYIII-YIl-lD/SII) 

AIN) =-.5 
8( 1 ) = S( 2 )/3. 

BI N ) = 1 . 

C( 1) = S( 2)/5. 

F( 1) = CYI 2 )-Y ( 1))/S(2)-Y1P 
FI N ) = 0. 

MI ) )=BU> 

SBI 1 )=Cl D/MU) 

G» 1 ) = FI D/WI 1) 

DO 30 1=2, N . 

MI 1 1 = B 1 1 )-At I )*SB I 1-1 ) 

SBI I)=C< 1 )/MI I) 

30 G< I) = IFI I )-A( I)*G( I-l) l/WII) 

FNI N ) =Gl N ) 

DO 40 1=2, N 
K=N+l-l 

40 EM I K )=GI < )-S.B(K)*EM(K+l) 

SLOPEl l)=-SI 2)/6.*(2.*EM(lHEM(2t)HYI2)-Y(l)t/Sm 
0050 1=2, N 

50 SIDPEl I)=SI l )/6.*(2.*EM(I >*EM( I-l ) ) ♦ I Y 1 1) -Yt 1-1 ) )/S II) 

IF I3.E3.14) WRITE 16,100) N, I XII) , YU ) , SLOPE (I) , EMII) , 1= l, N) 

100 FORMAT I2X15HN0. OF POINTS *I3/10X5HX 15X5HY 15XSHSL0PE15X 5H 
1EM /I4F20.8)) 

RETURN 

END 
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HBFTC St-.lNE DECK 

SUBR3U) INF SPLINE (X,Y,N, SLOPE , EM) 

DIMENSION X( 50) ,Y(50), SC 50) ,A(50) ,3 (50) ,C f 50) , F (50 ) ,W 150 ) , S B( 53 ) , 
1G(50),FM( 50) , SLOPF ( 50) 

COMMON Q 
INTEGER 0 
00 10 1 = 2, N 
10 S( I )=X( I)-X( I- 1) 

N0=N-1 

00 20 I=2*N3 
A( I)=S( I)/6. 

8( I )*( S( DtSII + 1) >/3. 

C( I)=S( I+l)/6. 

20 F( I )=( Y( I + l)-yf I) )/S( !♦!)-( YCI) -Y(l-l) )/SU) 

A( N) =-.5 
B( 1 )=I. 

B(N) = 1. 

Cf l)=-.5 
F( 11 = 0. 

F( N ) = 0» 

M 1 l=B« 1 » 

SB( 1)=C( 1 >/W( 1) 

G( l> = 0. 

DO 30 1=2, N 

M( I)=BII)-A( I }*Sb ( I — 1 1 
SB(I)-Cm/H(I) 

30 G( II = (F( T)-Am*GU-l) l/W(I ) 

EM(V)-C;\> 

DO AO 1=2, N 
K=N*1-I 

AO eVt;=Gm-SB(K)*EM(K+l) 

S'.OPL-m=-S( 2)/6.*(2.*EM(l )+6M(2))*(Y(2I-Ym )/S(2> 

0050 1=2, N 

50 SLOPF( I) = S(I ,/6.*(r.*EM(I )+EM( f-i) )♦ C YU I -Y( 1-1 ) ) /S ( 1 1 

IF (C.E3.13) WRITE ( 6 , 100 ) N, (XU) , VI 1 1 ,SLOPE(I> ,FHCn , I* I, N) 

100 FORMAT ( 2X15HNO. OF POINTS =I3/I0X5HX 15X5HY 15X5HSL0PE15X5H 
1EM / { AF20.8) ) 

RETURN 

END 




26 


SIBFTC INTGRL OECk 

SMRR3UHNF INTGRL (X,Y,N # SUM) 

DIMENSION X(50) ,Y<50) f S(50) ,A(50I ,9(50) f C <501 f F (53 ) t W (SO) ,S8f 501 , 

1 GC 50)rEH( 50 ) » SUM ( 50) 

COMMON SR W 
INTEGFR SR W 
00 1C I * 2 » N 

io s( u*xm-xc r-n 

N0=N-1 

DO 20 f*2*N3 

At I )“Sf n/6.0 

6( !)*< SC I)+S( 1 + 1) )/3.0 

cc n=st un/ 6 .o 

20 fc i ) = c yc i«n-Y(i ))/su + n-(Ycn-Yn-in/sm 

5 A<N)=-,5 

B( 11=1." 

BC N ) = 1 .0 

c <:>=-.5 

, F(1) = 0.0 

F(N) = 0.0 

I 

t. SB(n«cm/wm 

G( 1 1 = 0.0 
DO 30 l = 2*N 

i wc n*Bui-A<n*SB(i-n 

SBcn=ciM/wcn 

30 G(i: = (F(IHU)*GlM))/wm 
* EMC N) =G( N ) 

J no 40 I = 2« N 

| K=N+1- I 

l 40 EM(K)=Gm-SB(K)*EM(K + n 

$ SUM! 1 ) =0.0 

i DO 50 K = 2f N 

50 SUM ( < ) = 5UM(K-1H-S<K)*C YCKM-YCK-ll 1 /2. O-S(K) **3*< EM(K) ♦EM(K-l) 1/2 
14.0 

IFCSRW.E0.17) WRITF(<S f 1000) N f ( XCI) ♦ YCI) , SUMC I) , EMC 11 # I* 1 ,N) 



RETURN 

1000 FORMAT ( 17HK NO. OF POINTS =I3/10X5HX l5X5^ f Y 15X5HSUM 
1 13X10H2ND DERIV./C4F20. 8)1 

END 


I 

1 
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$IBFTC 1 ININT DECK 

SUBROUTINE L IN INT ( XI , Y1 , X , Y ,TN, MX, MY ,F ) 

COHHDN < 

DIMENSION X( MX ) ,Y(MY) ,TN(MX,MYI 
DO 10 J3 = l ,M X 

10 IFIXI.LE.XI J3HG0 TO 20 
J3=MX 

20 DO 30 J4=1,NY 
30 IFIY1 .LF.Y(J4)JG0 TO 40 
J4=MY 

40 J 1=J3-I 
J 2 *J 4- 1 

EPSl*IXl-X(jm/(XU3l-X(Jll> 

6PS2=(Yl-Y(J2n/I Yf J4I-YIJ2H 

EP$3=1.-EP$1 

FPS4=l.-EPS2 

F=TN< Jl, J2)*EPS3*EPS4+TN( J3,J2»*EPSl*EPS4*TNI Jl, J4)*EP$2*EPS3* 
ITN(J3,J4)*EP$1*EPS2 

IFK.EO. 14) WRITE (6»1) XI, Y1,F,JI,J?,FPS1 ,EPS2 
I FORMAT I3H LININT3F10.5,2I3,2F10.5» 

K=0 

RETURN 

END 


■j 





28 


SIBFTC CONTAL DECK 

SUBROUTINE CONTAL ( WA , WTFL » IND » I » WT 1 
DIMENSION SPEEO 1 3 1 , WEI GHT! 3) 

135 GO TO ( 140, 150,210,270, 370) ,!ND 
140 SPFFDI 1 1 * WA 
WEIGHT!!) = WTFL 
WA = WT/WTEL+WA 
TND * 2 
RETURN 

150 IF (IWTFL-WFIGHT1 1) ) /! Wft-SPEED <11 1 ) 180,180,160 

160 SP FED! 2 1 = WA 

WA = I WT-WTFL l/IWTFL- WEIGH Til) ) 

1 *( W A- SPEED! 1 ) 1+WA 

IF (ABSIWA—SPEEO (2)1-100.0) 166,166,161 

161 IFIWA-SPEED! 2)1163,163,162 

162 WA = SPEED! 2 1+100.0 
GO TO 166 

163 WA = SPFFO! ? 1- 100.0 
166 SPEED! 1) = SPEED! 2 ) 

WEIGHT! 11 = WTFL 
RFTURN 

170 WRITE (6,10001 I,WTFL 
IND = 6 
RETURN 
180 IND = 3 

IF (WTFL.GF.WT1 GO TO 140 
IF (SPEED! ll-WAl 190,200,200 
190 SPEED! 21 = SPEED! 11 

SPEED! 11 = ?.0*SPEEO( 1 l-WA 
SPEED! 31 = WA 
WEIGHT! 2) = WEIGHT! 1) 

WEIGHT! 31 = WTFL 
WA = SPEEDU1 
RETURN 

200 SPEED! 21 = WA 

SPEED! 31 = SPEED! 1 1 
SPEED! 11 = 2.0*WA-SPEFD(1) 

WEIGHT! 21 = WTFL 
WEIGHTO) = WEIGHT(l) 

WA « SPEED! 1 1 
RETURN 

210 WEIGHT! 11 = WTFL 

IF (WTFL.GE.WT1 GO TO 140 
IF ( WEIGHT! 1 l-WE IGHT( 2) 1 230,380,220 
220 WEIGHTO) = WEIGHT!?) 

WEIGHT! 2) * WEIGHTI11 
SPEED! 3) » SPEED! 2) 

SPEED! 2) = SPEED ( 1 ) 

' SPEED! 1) = 2.0* >PEE9 ( 2 1- SPEED I 31 
W'. * SPEED! 1) 

RETURN 

230 IF CSPFEDI31-SPEEDI1 1-10.0) 170,170,240 
240 IND = 4 

245 IF (WEIGHT! 31- WEIGHT 1 1) 1 260,260,250 
250 WA » ( SPEFOdl+SPEED! ?.) 1/2.0 
RETURN 
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260 WA = i SPFFDm+SPEEDl 2) )/2.0 
RETURN 

270 IF ( SPEED! 3)-SPEED!l J-10.0) 170,170,280 
280 IF ( WTFL -WEIGHT! 2) I 320,350,290 
290 IF (WA-SPEE012)) 310,300,300 
300 SPEED! 1) * SP EFD ! 2 ) 

SPEED! 2) * WA 
WFIGHTm = WFIGHT12I 
WEIGHT! 2) = WTFL 
GO TO 245 

310 SPFED! 3) = SPEED! 2 ) 

SPEED! 2 ) = WA 
WEIGHT! 3) * WE IGHT! 2) 

WEIGHT! 2) = WTFL 
GO TO 245 

320 IF IWA-SPEED12I) 340,330,330 
330 WEIGHT! 3) = WTFL 
SPEED! 3) = WA 
GO TO 245 

340 WEIGHT! 1) = WTFL 
SPEED! 1) = WA 
GO TO 245 
350 IND = 5 

IF ! WA-SPEED 121) 380,360,360 
360 SPEED! 1) = SPEED! 2) 

WEIGHT! 1 ) = WE IGHT! 2 ) 

SP EED! 2; = ! SPEED ! 1 ) +SPEED 1 3) ) / 2. 0 
W A = SPEED! 21 
RETURN 
370 IND = 4 

WEIGHT! 2) = WTFL 
WA = ! SP EFD! I ) +SPFE0! 2)1 /2, 0 
RETURN 
380 IND * 5 

390 WEIGHT! 3) * WEIGHT(2) 

SPEED! 3) = SPEED! 21 

SP FED! 2) = (SPEED ! 1 1 + SPEED ( 3) ) /?• 

WA * SPEED! 2 ) 

RETURN 

1000 FORMAT (/12H FIXED LINE I2,12H, MAX WT * F10.6I 
END 
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SAMPLE OUTFIT 

For aid in checking the operation of the program and as an illustration 
of the results obtained by the program, sample output is included below. 

The case used here is the same radial gas turbine rotor with splitter blades 
that was used for the numerical example in reference (2). The first out- 
put is from the calling program (Q3D) and consists of the values of variables 
MOB, H, M, ana S. The value of MOB determines whether the meridional solu- 
tion (FIXED) or the blade-to-blade solution (3PLBLD) or both is desired. In 
the example below MOB = 2 which gives both solutions. The value of H, M, r 
S is used to indicate whether the hub, meridional or shroud blade-to-blade 
solution is desiredj the value 1 indicates that particular solution is desired, 
the value 0 not. In the example H - 1 , and M and S = 0. 

The data given immediately following "RUN NO. 1 INPUT DATA CARD LISTING", 
up to the line before "STAG. SPEED OF SOUND AT INLET", are a listing of the 
data on the input data cards for the meridional solution. This is followed 
by a list of the maximum computed streamline change at each iteration for the 
meridional solution. Convergence to the specified accuracy of .001" was 
obtained after 105 iterations. Then the calculated information is printed 
out for every 5th streamline. Data at more streamlines will be printed out 
if the value of NPRT is reduced from 5. 

The output for the meridional solution is given in two parts. The first 
part gives values of DN, Z, R, WA, PRESS, WTR, and RC and the maximum stream- 
line change for iteration No. IOC. The second part for iteration No. 107 
gives values for ALPHA, RC, 31, BETA, TT, SA, SB, SC, and SD, followed by a 
repetition of the values of DN, Z, R, WA, PRESS, WTR, and RC. The next out- 
put U the "STREAMLINE SPACING ALONG NORMAL", "BLADE COORDINATES", and "IN- 
LET ANGLES", which is computed by the meridional program (FIXED) and is used 
for input to the blade-to-blade program (SPLELD). 

Next is the output from the blade-to-blade program at the hub. The first 
line is "RUN NO. 2 INPUT DATA CARD LISTING". The next two lines are data f'om 
input cards. The following lines up to the 3rd line before "THETa- CALCULATED 
AND/OR II °UT", are computed by the program, and are the data which would have 
to be supplied if not furnished by the meridional program. The second line 
before "THETA- CALCULATED AND/OR INPUT" is data from an input card with a 1 
in card column 5, indicating further data is given. A blank card here would 
indicate no more input for this blade-to-blade solution. The line just 
before "THETA- CALCULATED AND/OR INPUT" is data from ai: input card giving 
information of splitter blade coordinates near the trailing edge that were 
altered from those computed by the meridional program. The next data are 
the Q coordinates of the Initial equally spaced streamlines at every 
quasi -orthcngonal. This is followed by the stagnation speed of sound, and 
then a list of the maximum predicted streamline change at each iteration 
lcr the blade-to blade solution. Convergence to the desired accuracy of 
.001" was obtained in 27 iterations. The desired output is printed out for 
every 5th streamline at each quasi-orthogonal. Again, data at more stream- 
lines can be obtained with a different value of NPRT. 

The mean and shroud blade-to-blade solutions would have also been obtained 
on the same computer run if we had set M and S equal to 1, and supplied the 
appropriate input cards. 
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INPUT DATA CARO LISTING 

4030.0000 0.6110 11.0000 

1950.0000 22 1.5000 0.02W 

2.0000 l.OOOC -25.1000 

1.1200 
2.5500 
0.3170 
2.5500 
2.2420 
2 . 1 * T 0 
1.4700 
0.7500 
- 0 . 


PUN NO . 1 


14 21 

7 7 

0 2 

0 11 

19 5 

2 0 

0.7590 

C. 8100 

1.5500 

1.8030 

0 . 

-0.0030 

1.3650 

1 • 7650 

3.0100 

2.7360 

2.1470 

2.1470 

3.0100 

2.6000 

0.7970 

0. 752t 

-0. 

-0. 

-0.0012 

-0.0044 

-0.6231 

-0.8604 

-0.2000 

0. 

1.1000 

1.2000 

2.4000 

2.600 > 

0.1860 

I 710 

0.1660 

0. 15* 0 

0.1460 

0. 1310 

0.1260 

c.nio 

0.1050 

0.0910 

0.0855 

0.0710 

0 3670 

0.0510 

-0.0040 

0. 4000 

0.6000 

1.0000 


0.8900 

1.0150 

2 1000 

2* 3500 

0.0300 

0.1400 

2m 1000 

2.3500 

2.5130 

2.32 30 

2.1470 

2. 1470 

2.2300 

1.8220 

0.7500 

0. 7500 

0 . 

-Oo 

0.0201 

-0. 0567 

1.0978 

-1.335? 

0.2000 

0.4C00 

1.4000 

1.6000 

2.8000 

3.0000 

0.1560 

0.1410 

0.1360 

0.1210 

o.: 160 

0.1010 

G. 0960 

0.0820 

0.0760 

0.0620 

0.0560 

0. 0420 

0. 3360 

0.0220 

0.8000 

1.2000 

1.4000 

1. 8000 


-0.1248 
-1.5725 
0.60D0 
1.8000 

3.2000 
0.1270 0.1160 

0.1065 0.0960 

0.0870 0.0760 

0.0670 0.0560 

0.0480 0.3360 

0.0280 0.0163 

O.OOPO -0..0050 

1.6000 2.0003 

2.2000 2.6333 


0.000 

1.3520 
3. 1530 
0.9430 
3.1500 
2.1600 
2.1470 
0.9230 
0.7500 
- 0 . 

-0.3966 

1.0300 

2.2000 

0.0330 

0.0300 

0.0270 

0.0240 

0.0220 

0.0190 

0.0160 

2.403a 

3.0200 


1.6667 1245.0300 
0.3313 0.2500 

0.0500 
1.2300 
2.9503 
0.5903 
2.3503 
2.1383 
2.1473 
1.1583 
0.7503 
- 0 . 

-0 >2339 

0.3303 
2.0300 


STAG. SPEED OF 

SOUND AT 

INLET 

■ 2011.55 




.2073 

ITERATION 

NO. 

1 

MAX. 

STREAMLINE 

CHANGE 

3 

0.397387 

ITERATION 

NO. 

2 

MAX. 

STREAMLINE 

CHANGE 

3 

0.364558 

ITERATION 

NO « 

3 

MAX. 

STREAMLINE 

CHANGE 

3 

0.338418 

ITERATION 

NO. 

4 

MAX. 

STREAMLINE 

CHANGE 

3 

0.317072 

ITERATION 

N3 . 

5 

MAX. 

STREAMLINE 

CHANGE 

3 

0.294279 

ITERAT 'ON 

NO. 

6 

MAX. 

STREAMLINE 

CHANGE 

3 

0.277526 

ITERATION 

NO. 

7 

MAX. 

STREAMLINE 

CHANGE 

a 

0.261206 

ITERATION 

NO. 

8 

MAX. 

STREAMLINE 

CHANGE 

3 

0.245780 

ITERATION 

NO. 

9 

M4X. 

STREAMLINE 

CHANGE 

3 

0.^31244 

ITERATION 

NO. 

10 

MAX. 

STREAMLINE 

CHANGE 

3 

0.217489 

ITERATION 

NO. 

11 

MAX. 

streamline 

CHANGE 

3 

0.204468 

ITERATION 

NO. 

12 

MAX. 

STREAMLINE 

CHANGF 

3 

0.192163 

ITERATION 

NO. 

13 

MAX. 

STREAMLINE 

CHANGE 

3 

0.110295 

ITERATION 

NO. 

14 

MAX. 

streamline 

CHANGE 

3 

0.169344 

ITERATION 

NO. 

15 

MAX. 

STREAMLINE 

CHANGE 

3 

3.159114 

ITERATION 

NO. 

16 

MAX. 

STREAMLINE 

CHANGE 

S 

0.149380 

ITERATION 

NO. 

17 

MAX. 

STREAMLINE 

CHANGE 

3 

0.140383 

IfCRATION 

NO * 

18 

MAX. 

streamline 

CHANGE 

- 

0.132053 

ITERATION 

NO. 

19 

MAX. 

streamline 

CHANGE 

3 

0.124201 

ITLPATION 

NO. 

?0 

MAX. 

streamline 

CHANGE 

3 

0.115800 

ITERATION 

NO. 

21 

MAX. 

streamline 

CHANGE 

3 

0.109090 

ITERATION 

NO. 

22 

MAX. 

STREAMLINE 

CHANGE 

3 

0.099879 

ITERATION 

NO. 

23 

MAX. 

STREAMLINE 

CHANGE 

3 

0.093197 

ITERATION 

NO. 

24 

MAX. 

STREAMLINE 

CHANGE 

3 

0.05751 1 

ITERATION 

NO. 

?5 

MAX. 

STREAMLINE 

CHANGE 

3 

0.032277 

ITERATION 

NO. 

26 

MAX. 

STREAMLINE 

CHANGE 

3 

0.077409 

ITER A! ION 

NO. 

27 

MAX. 

STREAMLINE 

CHANGE 

« 

0.072972 

ITERATION 

NO. 

28 

MAX. 

STREAMLINE 

CHANGE 

3 

0.068739 

ITERATION 

NO. 

29 

*AX. 

STREAMLINE 

CHANGE 

3 

0.064801 

ITERATION 

NO. 

30 

MAX. 

STREAMLINE 

CHANGE 

3 

0.061076 
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ITERATION 

NO • 

01 

MAX. 

STREAMLINE 

CHANGE 

s 

0*057551 

ITERATION 

NO. 

02 

MAX. 

STREAMLINE 

CHANGE 

s 

0.054177 

ITERAt ION 

NO. 

33 

MAX. 

STREAMLINE 

CHANGE 

s 

0.051032 

ITERATION 

NO. 

34 

MAX. 

STREAMLINE 

CHANGE 

= 

0.043114 

ITERATION 

NO. 

35 

MAX. 

STREAMLINE 

CHANGE 

s 

0.045275 

ITERATION 

NO. 

36 

MAX. 

streamline 

CHANGE 


5.042680 

ITERATION 

NO. 

37 

MAX. 

streamline 

CHANGE 


0.040203 

ITERATION 

NO. 

38 

MAX. 

STREAMLINE 

CHANGE 


0.037993 

ITERATION 

NO. 

39 

MAX. 

STREAMLINE 

CHANGE 


0.035860 

ITERATION 

NO. 

40 

MAX. 

STREAMLINE 

CHANGE 


0.033894 

ITERATION 

NO. 

41 

MAX. 

STREAMLINE 

CHANGE 

= 

0.031966 

ITERATION 

N'l • 

42 

MAX. 

streamline 

CHANGE 

sr 

0.030159 

ITERATION 

NO. 

43 

MAX. 

STREAMLINE 

CHANGE 

3 

0.028456 

ITERATION 

NO. 

44 

MAX. 

STREAMLINE 

CHANGE 

= 

0.026857 

ITERATION 

NO. 

45 

MAX. 

STREAMLINE 

CHANGE 

= 

0.025363 

ITERATION 

NO. 

46 

MAX. 

STREAMLINE 

CHANGE 

- 

0.023967 

ITERATION 

NO. 

47 

MAX. 

STREAMLINE 

CHANGE 

= 

0.022629 

ITERATION 

NO. 

48 

MAX. 

STREAMLINE 

CHANGE 

3 

0.021369 

itera; ION 

NO. 

49 

MAX. 

STREAMLINE 

CHANGE 

= 

0.020179 

ITERATION 

NO. 

50 

MAX. 

STREAMLINE 

CHANGE 

3 

0.019052 

ITERATION 

NO. 

51 

MAX. 

STREAMLINE 

CHANGE 

3 

0.017998 

ITERATION 

NO. 

52 

MAX. 

STREAMLINE 

CHANGE 

3 

0.017001 

ITERATION 

NO. 

53 

MAX. 

STREAMLINE 

CHANGE 

3 

0.016058 

ITERATION 

NO. 

54 

MAX. 

STREAMLINE 

CHANGE 

3 

0.015168 

ITERATION 

NO. 

55 

MAX. 

STREAMLINE 

CHANGE 

3 

0.014330 

ITERATION 

NO. 

56 

MAX. 

streamline 

CHANGE 

S 

0.013539 

ITERATION 

NO. 

57 

MAX* 

STREAMLI NE 

CHANGE 


0.012791 

ITERATION 

NO. 

58 

MAX. 

STREAMLINE 
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3 

0.012087 

ITERATION 

NO. 

59 

MAX. 

STREAMLINE 

CHANGE 

3 

0.011487 

ITERATION 

NO. 

60 
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STREAMLINE 

CHANGE 

3v 

5.010813 

ITERATION 

NO. 

61 

MAX. 

STREAMLINE 

CHANGE 

3 

0.010206 

ITERATION 

NO. 

62 
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STREAMLINE 
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3 

0.009696 
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NO. 

63 
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STREAMLINE 
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= 

0.009128 

ITERATION 

NO. 

64 
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STREAMLINE 
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3 

0.008669 
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NO. 

65 
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3 

0.008158 

ITERATION 
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66 
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STREAMLINE 
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3 
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ITERATION 

NO. 

67 
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S.REAMLINE 
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S 

0.007291 

ITERATION 

*3. 

58 
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STREAMLI NE 
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3 

0.006937 

ITERATION 

NO * 

69 

MAX. 

streamline 
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3 

0.006519 

ITERATION 

NO. 

70 
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3 

0.006204 

ITERATION 

NO. 

71 
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3 

0.005890 

ITERATION 

NO. 

72 
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STREAMLINE 
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S 

0.005517 

ITERATION 

NO. 

73 
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STREAMLINE 
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3 

0.005248 

ITERATION 
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74 
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STREAMLINE 
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3 

0.004985 

ITERATION 

NO. 

75 
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STREAMLINE 
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S 

0.004665 
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NO. 

76 
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STREAMLINE 
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3 

0.004424 

ITERATION 

NO. 

77 
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STREAMLINE 
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3 

0.004201 
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78 

MAX. 

streamline 
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3 

0.003997 

ITERATION 

NO. 

79 
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STREAMLINE 
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3 

0.033801 
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NO. 

80 

MAX. 

STREAMLINE 
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3 

0.003538 

ITERATION 

NO. 

81 

MAX. 

STREAMLINE 

CHANGE 

S 

0.003360 

ITERATION 

NO. 

82 

MAX. 

STREAMLINE 

CHANGE 

3 

0.003189 

ITERATION 

NO. 

a:. 

MAX* 

STREAMLINE 

CHANGE 

3 

0.003028 

ITERATION 

NO. 

84 

MAX. 

STREAMLINE 

CHANGE 

S 

0 .002872 

ITERATION 

NO. 

85 

MAX. 

STREAMLINE 

CHANGE 

3 

0.002728 

ITERATION 

NO. 

86 

MAX. 

STREAMLINE 

CHANGE 

9 

0.002603 

ITERATION 

NO. 

87 

MAX. 

STREAMLINE 

CHANGE 

3 

0.002406 

ITERATION 

NO. 

88 

MAX. 

STREAMLINE 

change 

9 

0.002284 

ITERATION 

NO. 

89 

MAX. 

STREAMLINE 

CHANGE 

3 

0.002178 

ITERATION 

NO. 

90 

MAX. 

STREAMLINE 

CHANGE 

3 

0.002058 
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RUN NO. 2 INPUT OATA CARO LISTING 


2? 8 

0 . 

16 22 

n 2 

3.3100 
8 0 

0.1500 

0.0330 

4030.0000 

0.0010 
0.03 05 

0.2003 

11.0000 

1.6567 

1245.0000 

? ? 

" 1 2 

1950.0300 

221.5000 

0.0247 

0.3310 

0.2500 6941.2549 

0 5 

-0.5661 

2 0 
-0.5600 

-24.5566 

-0.5548 

G. 0000 
-0.5482 

0.1500 

-0.5384 

-0.5268 

-0.5125 

-0.4984 

-0.5103 

-0.6036 

-0.80 85 

-1.0977 

-1.3724 

-1.7284 

-2.0844 

0.0151 

-0,0132 

-0.0164 

-0.0229 

-0.0328 

-0.0444 

-0 .0587 

-0.0729 

-0.0972 

-0. 1918 

-0.3765 

-0.6005 

-0.8012 

-U1572 

-1.513? 

-0.2705 

-0.2968 

-0.3020 

-0.3005 

-C.3184 

-0.3300 

-0.3443 

-0.3585 

-0.2705 

-0.2744 

-0.2692 

-0.2626 

-0*2528 

-0.2412 

-0.2269 

-0.2128 

0. 

0 . 

-0.0002 

0.0025 

0.0117 

0.0264 

0,049? 

0.0786 

0. U?9 

0.1471 

0.1750 

0.1958 

0.2125 

0.2375 

0.2625 

0 • 2758 

0.2508 

0.2167 

0.1858 

0.1518 

0.1225 

0.0965 

0.0769 

0.0664 

0.0627 

0.0625 

0.0625 

0.0625 

0.0625 

0.0625 

v.0032 

0.0034 

0.0043 

0.0057 

0.0079 

0.0109 

0,0142 

0.0172 

0.0167 

0.0146 

0.0112 

0.0092 

0.0085 

0.0083 

0.0082 


l _ . „ 

711-0* 33 890 311-0.31900 712-0.23630 012-0*28660 812-0.286600 
THFT A-C ALCULATEO ANO/OR INPUT 


-0.5561 

-0.5600 

-0.5548 

-0.5482 

-0.5384 

-0.5268 

-0.5125 

-0.4984 

-2.0844 

-0.5103 

-0.6036 

-0. 8085 

-1.0977 

- 1.3T24 

■1.7294 

-0.5275 

-0.5337 

-0.5295 

-0.5243 

-0.5164 

-0.6071 

-0.4957 

-0.4844 

-2.0558 

-0.4897 

-0.5830 

-0.7869 

-1.0728 

-1.3438 

-1.6998 

-0.4990 

-0.5074 

-0.5042 

-0.5003 

-0.4944 

-0.4974 

-0.4789 

-0.4704 

-2.0273 

-0.4690 

-0.5624 

-0. 7653 

-1.0479 

-1.3X53 

-1.6713 

-0.4704 

-0.4810 

-0.4789 

-0.4763 

-0.4724 

-0,4678 

-0.4621 

-0.4564 

-1.9987 

-0.4484 

-0.5419 

-0.7437 

-1.0231 

-1.2867 

-1.6427 

-0.4418 

-0.4547 

-0,4537 

-0.4524 

-0.4504 

-0.4481 

-0.4453 

-0.4424 

-1.9701 

-0.4277 

-0.5213 

-0.7221 

-0.9982 

-1.2582 

-1.6141 

-0.4133 

-0.4284 

-0.4284 

-0.42 84 

-0.4284 

-0.4284 

-0.4204 

-0.4284 

-1.9416 

-0.4071 

-0.5007 

-0. 7005 

-0.973 V 

-1.2296 

-1.5856 

-0.3847 

-0.4021 

-0.4031 

-0.4044 

-0.406V 

-0.4087 

-0.4116 

-0.4145 

-1.9130 

-0.3864 

-0.480* 

-0.6789 

-0. 9485 

-1.2013 

-1.5570 

-0.3562 

-0.3757 

-".3779 

-0.3805 

-0.384V 

-0.3890 

-0.3948 

-0.4005 

-1.8845 

-0.3657 

0.4595 

-0.6573 

-0.9236 

-1.1725 

-1,5285 

-0.3276 

-0.3A9* 

-0.3526 

-0.3565 

-0.3624 

-0.3694 

-0.^780 

-0.3865 

-1.8559 

-0. 3451 

- 0. 4 £ 89 

-0.6357 

-0.8988 

-1.1439 

-1.4999 

-0.2990 

-0.3231 

-0.3273 

-0.3325 

-0.3404 

-0.3497 

-0.3612 

-0.3725 

-1.8273 

-0.3244 

-0.41^3 

-0.6141 

-0.8739 

-1.1154 

-1.4713 

-0.2705 

-0.2968 

-0.3020 

-0.3085 

-0.3184 

-0.3300 

-0.3389 

-0.3190 

-1.7988 

-0.3038 

-0.3977 

-0.5925 

-0.CV91 

-1.0869 

-1.4428 



40 


-U.2705 

-0.2744 

-0.2692 

-0.2626 

-0.2528 

-0.2412 

-0.2363 

•0.2866 

•1.7988 

-0.3038 

-0.3977 

-0.5925 

-0.8491 

-1 .0868 

-1.4428 

•0.2419 

-0.2481 

-0.2439 

-0.2387 

-0.2308 

-0.2215 

-0.2101 

-0.1988 

-1.7702 

-0.2331 

-0.3771 

-0.5709 

•0.8247 

-1.0582 

-1.4142 

-0.2134 

-0.2218 

-0.2186 

-0.2147 

•0.2088 

-0.2013 

-0.1933 

-0.1848 

-1.7417 

-0.2625 

-0.3565 

-0.5453 

-0.7993 

-1.0297 

-1.3857 

-0.1848 

-0.1954 

-0.193.* 

-0.1907 

-0.1858 

-0.1822 

-0.1765 

-0.1708 

-1.7131 

-0.2418 

-0.3359 

-0.5277 

-0.7745 

'1.0011 

-1.3571 

-0.1562 

-0.1691 

-0.1681 

-0.1668 

-0.1648 

-0.1625 

-0.1597 

-0,1568 

-1.6845 

-0.2212 

-0.3153 

-0.5061 

•0.7496 

-0.9726 

-J.3285 

-0.1277 

-0.1428 

-0.142G 

-0.1426 

-0.1428 

-0.1428 

-0.142P 

-0.1428 

-1.6660 

-0.2005 

-0.2947 

-0.4 845 

-0.7248 

-0.9440 

-1.3000 

-0.0991 

-0.1165 

-0.1175 

-0.1188 

-0.1208 

-0.12^1 

-0.1260 

-0.1289 

-1.6274 

-0.1799 

-0.2742 

-0. 4629 

-0.5999 

-0.9154 

-1.2714 

-0.0706 

-0.0901 

-0.0923 

-0.0949 

-0.0988 

-0.1034 

-0.1092 

-0.1149 

-1.5989 

-0.1592 

-0.2536 

-0.4413 

-0.6750 

-0.8369 

-1.2429 

-0.0*20 

-0.0638 

-0.0670 

-0.0709 

-0.0768 

-0.0938 

-0.0924 

-0.1009 

-1.5703 

-0.1386 

-0.2330 

-0.4197 

-0.6502 

-0.8533 

-1.2143 

-0.0134 

-0.0375 

-0.0417 

-0. 0469 

-0.0548 

-0,0641 

-0.0756 

-0.0869 

-1.5417 

-0.1179 

-0.2124 

-0.3981 

-0.8253 

-0.9298 

-1.1857 

0.0151 

-0.011? 

-0.0164 

-0.0229 

-0.0328 

-0.0444 

-0.0587 

-0.0729 

-1.5132 

-0.0972 

-0.1918 

-0.3765 

-0.6005 

-0.8012 

-1.1572 


STAG, speed of 

SOUND AT 

INLET 

* 2011.55 




ITERATION 

NO. 

1 

MAX. 

STREAMLINE 

CHANGE 

JS 

0. 068 444 

ITERATION 

NO. 

2 

MAX. 

STREAMLINE 

CHANGE 

3 

0.058218 

ITERATION 

NO. 

3 

MAX. 

STREAMLINE 

CHANGE 

3 

0.049416 

ITERATION 

NO. 

4 

MAX. 

STRF AMLINE 

CHANGE 

- 

0.041949 

ITERATION 

NO. 

5 

MAX. 

STREAMLINE 

CHANGE 

3 

0.035614 

ITERATION 

NO. 

6 

MAX. 

STREAMLINE 

CHANGE 

a 

0.030236 

ITERATION 

NO. 

7 

MAX. 

STREAMLINE 

CHANGE 

- 

0.025670 

ITERATION 

NO, 

8 

MAX. 

STREAMLINE 

CHANGE 

3 

0.021794 

ITERATION 

NO. 

9 

MAX. 

STREAMLINE 

CHANGE 

3 

0.018503 

ITERATION 

NO. 

10 

MAX. 

STREAMuINE 

CHANGE 

3 

0.015709 

ITERATION 

NO. 

11 

MAX* 

STREAMLINE 

CHANGE 

* 

0.013338 

ITERATION 

NO. 

12 

MAX. 

STREAMLINE 

CHANGE 

3 

0.011324 

ITERATION 

NO. 

1 3 

MAX. 

STREAMLINE 

CHANGE 

3 

0.009614 

ITERATION 

NO. 

14 

MAX. 

STREAMLINE 

CHANGE 

3 

0.008163 

ITERATION 

NO. 

15 

MAX. 

STREAMLINE 

CHANGE 

3 

0.006931 

ITERATION 

NO. 

16 

MAX. 

STREAMLINE 

CHANGE 

* 

0.005885 

ITERATION 

NO. 

17 

MAX. 

STREAMLINE 

CHANGE 

3 

0.004996 

ITERATION 

NO. 

18 

MAX. 

SUEAMLINE 

CHANGE 

3 

0.004248 

ITERATION 

NO. 

19 

MAX. 

STRF AMLI NE 

CHANGE 

3 

0.003601 

ITERATION 

NO. 

?0 

MAX. 

STREAMLINE 

CHANGE 

3 

0.003061 

ITERATION 

NO. 

21 

MAX. 

streamline 

CHANGE 

3 

0.002595 

ITERATION 

NO. 

22 

MAX. 

STREAMLINE 

CHANGE 

3 

0.002207 

ITERATION 

NO. 

23 

MAX. 

STREAMLINE 

CHANGE 

3 

0.001876 

ITERATION 

NO. 

24 

MAX. 

STREAMLINE 

CHANGE 

3 

0.001588 



ITERATION NO. 25 MAX. STREAMLINE CHANGE * 3.001350 

ITERATION NT. 26 HA X* STREAMLINE CHANSE * 3.031148 

ITERATION NT. 27 MAX. STREAMLINE CHANGE * 3.900976 

ITERATION NT. 28 MAX. STREAMLINE CHANSE » 3.000829 
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